DUDLEY  KNOX  LIBRARY 
NAVAL  POSTGRADUATE  SCHOOt 
SqnTEREY  CA  93943-5101 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OMB  No    0704-0188 


Public  reporting  burden  tor  this  collection  of  information  is  estimated  to  average  I  hour  per  response,  including  the  time  tor  reviewing  instructions,  searching  existing  data  sources, 
gathering  and  maintaining  the  data  needed,  and  completing  and  reviewmq  the  collection  of  information  Send  comments  regarding  this  burden  estimate  or  anv  other  aspect  of  this 
collection  of  information,  including  suggestions  for  reducing  this  burden  to  Washington  Headauarters  Services.  Directorate  tor  information  Operations  and  Reports.  1215  Jefferson 
Davis  Highway.  Suite  1204.  Arlington.  VA   22202-4302.  and  to  the  Office  of  Management  and  Budget.  Paperwork  Reduction  Proiect  (0704-0188).  Washington.  DC  20503 


1.    AGENCY  USE  ONLY  (Leave  blank) 


2.  REPORT  DATE 

March  1993 


3.  REPORT  TYPE    AND  DATES  COVERED 
Master's  Thesis     March  1993 


4.  TITLE  AND  SUBTITLE 


Theoretical    Study  of  Laminar    Film  Condensation    on  Horizontal 
Elliptical    Tubes  Under  Conditions    of  Free  and  Forced  Convection 


6.  AUTHOR(S) 

Vance  Hiro  Adams 


5.    FUNDING  NUMBERS 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Naval  Postgraduate    School 
Monterey,  CA    93943-5000 


8.    PERFORMING  ORGANIZATION 
REPORT  NUMBER 


9.  SPONSORING  /MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

Naval  Postgraduate    School 
Monterey,  CA    93943-5000 


10.  SPONSORING/MONITORING 
AGENCY  REPORT  NUMBER 


11.  SUPPLEMENTARY  NOTES 

The  views  expressed   are  those  of  the  author   and  do  not  reflect   the  official    policy  or  position 
of  the  Department    of  Defense   or  the  U.S.  Government 


12a.  DISTRIBUTION /AVAILABILITY  STATEMENT 

Approved   for  public  release:   distribution    is  unlimited 


12b.  DISTRIBUTION  CODE 


13.  ABSTRACT  (Maximum  200  words) 

Analytical   studies  have  been  made  of  laminar  film  condensation    on  a  horizontal   elliptical 
cylinder  in  a  pure  saturated    vapor  under  conditions  of  free  and  forced  convection.    Estimation 
of  interfacial    shear  stress  was  made  in  two  ways:  the  first  involving   an  asymptotic   value  of  the 
shear    stress    under    conditions     of    infinite     condensation      rate    and    the    second    based    on 
simultaneously    solving   the  two-phase  vapor   boundary    layer   and  condensate    equations.     The 
latter  approach   enables  the  determination    of  the  vapor  boundary    layer  separation    point.    For 
the  assumption   of  asymptotic    shear  stress,  effects  of  surface   tension  and  pressure  gradient   in 
the  condensate    film  have  been  included.     At  the  extremes  of  eccentricity,    corresponding    to  a 
circular   tube  and  a  vertical   plate,  the  results  are  compared   with  theoretical    and  experimental 
work  of  others.     Improvement    in  the  condensation    heat  transfer    coefficient    was  found    for 
elliptical    tubes  under  both  free  and  forced  convection    conditions    when  compared   to  circular 
tubes  of  the  same  surface    area.     In  the  latter   case,  this  improvement    was  due  mainly   to  the 
reduced  drag  of  the  elliptical   tube  providing  a  higher  vapor  velocity  for  the  same  pressure  drop 
as  that  across  a  circular   tube. 


14.  SUBJECT  TERMS 

Laminar   Film  Condensation,    Horizontal    Elliptical    Tube, 
Free  and  Forced  Convection 


15.  NUMBER  OF  PAGES 

124 


16.  PRICE  CODE 


17.    SECURITY  CLASSIFICATION 
OF  REPORT 

unclassified 


18.    SECURITY  CLASSIFICATION 
OF  THIS  PAGE 

unclassified 


19 


SECURITY  CLASSIFICATION 
OF  ABSTRACT 

unclassified 


20.  LIMITATION  OF  ABSTRACT 

SAR 


NSN  7540-01-280-5500 


Itand 


14% 


Approved   for  public  release,   distribution    unlimited 

THEORETICAL    STUDY   OF  LAMINAR  FILM  CONDENSATION 

ON  HORIZONTAL    ELLIPTICAL   TUBES  UNDER  CONDITIONS 

OF  FREE  AND  FORCED   CONVECTION 


by 


Vance  Hiro  Adams 

Lieutenant    Commander,    United   States   Navy 

B.S.,  United  States   Naval  Academy,   1978 

Submitted    in  partial    fulfillment    of  the 
requirements    for  the  degrees  of 


MASTER  OF  SCIENCE  IN  MECHANICAL  ENGINEERING 

and  the  degree  of 
MECHANICAL  ENGINEER 

from  the 

NAVAL  POSTGRADUATE    SCHOOL 

March  1993 


ABSTRACT 

Analytical  studies  have  been  made  of  laminar  film  condensation  on  a  horizontal 
elliptical  cylinder  in  a  pure  saturated  vapor  under  conditions  of  free  and  forced  convection. 
Estimation  of  interfacial  shear  stress  was  made  in  two  ways:  the  first  involving  an  asymptotic 
value  of  the  shear  stress  under  conditions  of  infinite  condensation  rate  and  the  second  based 
on  simultaneously  solving  the  two- phase  vapor  boundary  layer  and  condensate  equations.  The 
latter  approach  enables  the  determination  of  the  vapor  boundary  layer  separation  point.  For 
the  assumption  of  asymptotic  shear  stress,  effects  of  surface  tension  and  pressure  gradient  in 
the  condensate  film  have  also  been  included.  At  the  extremes  of  eccentricity,  corresponding 
to  a  circular  tube  and  a  vertical  flat  plate,  the  results  are  compared  with  theoretical  and 
experimental  work  of  others.  Improvement  in  the  condensation  heat  transfer  coefficient  was 
found  for  elliptical  tubes  under  both  free  and  forced  convection  conditions  when  compared 
to  circular  tubes  of  the  same  surface  area.  In  the  latter  case,  this  improvement  was  due  mainly 
to  the  reduced  drag  of  the  elliptical  tube  providing  a  higher  vapor  velocity  for  the  same 
pressure   drop  as  that  across  a  circular    tube. 
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NOMENCLATURE 

a  semi-major  axis  dimension 

Ac  tube  cross-sectional   area 

b  semi-minor  axis  dimension 

B:  as  defined    by  Equation    (11-35) 

B-,  as  defined   by  Equation    (11-52) 

Bo  Bond  number  (p(gDe-  /  o) 
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g  gravitational     constant 
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hf  latent   heat  of  evaporation 

k  elliptical    tube  eccentricity    (b/a) 

m  condensate    mass  flux  rate 

Nu  Nusselt  number  (ocDe/A.,) 

p  fluid  pressure 

p^  fluid  pressure   due  to  potential    flow 

p0  fluid  pressure   due  to  surface   tension 

P  dimensionless    parameter    defined   by  Equation    (11-85) 

r  radial   dimension,   cylindical    coordinate    system 
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R(<j>)  radius   of  curvature 

Re  two-phase  Reynolds    number   (p^^/ r^) 

Rev  vapor   Reynolds    number  (pvU00De/r|v) 

Twall  tube  wall  temperature 

Tsat  vapor  saturation    temperature 
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U  vapor    boundary- layer    streamwise    velocity 
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U^  dimensionless    vapor  potential    velocity   defined   by  Equation    (11-108) 

U.,  vapor  free-stream   velocity 

v  condensate    film  normal   velocity 
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I.    INTRODUCTION 

A.        MOTIVATION    FOR  INVESTIGATION     OF  ELLIPTICAL    TUBES 

For  many  years,  condenser  design  in  the  United  States  has  been  based  on  the  Heat 
Exchange  Institute  (H.E.I.)  or  Tubular  Exchange  Manufacturers  Association  (T.E.M.A.) 
standards  [1,2]  which  rely  on  tube-side  empirical  relations  and  averaged  data  to  predict 
heat  transfer  parameters.  Since  there  is  a  large  degree  of  uncertainty  in  the  accuracy 
of  the  prediction,  the  condensers  are  over  designed  to  ensure  sufficient  margin  of 
reliability.  The  result  is  an  excessively  large,  expensive  condenser  for  the  desired 
thermal  duty  rating  requirements.  The  majority  of  condensers  currently  used  in  the 
Navy  were  designed   under  these  standards. 

Within  the  last  two  decades,  computer  modeling  has  been  used  to  more  accurately 
predict  the  heat- transfer   coefficients    of  condensers.    There  are  two  basic  approaches   to 
computer  modelling,  an  integrated    and  a  differential    approach,  of  which  the  former  is 
most  commonly  used  in  condenser  design  (Walker  [3]).  Integrated   approaches   divide  the 
condenser  into  zones  with  associated   mean  heat  transfer   and  fluid  properties  which  are 
typically   determined   from  semi-empirical  relations.    These  zones  are  integrated   over  the 
whole  condenser  and  provide  overall   performance    predictions    with  respect  to  shell  side 
pressure   drop  and  duty  rating.    Though   easy  to  use,  integrated    methods  are  restrictive 
when    used    in    conditions     which    are    outside    of   the    fluid    parameter     and    geometry 
constraints     of    the    empirical     relations     employed     in    the    approach.        Differential 
approaches   solve  the  basic  equations   of  fluid  mechanics   and  heat  transfer  and  provide 
an  understanding    of  local  variations    within  the  condenser,  but  are  too  complex  for  use 
in  condenser    design.     The   improvements    in  accuracy    as  a  rsult   of  computer    design 


methods  are  a  function  of  the  empirical  or  theoretical    models  used  in  their  development. 
Marto  [4]  provides   a  comprehensive    overview   of  the  state  of  computer   design  methods. 

The  purpose  of  empirical  and  theoretical  research  into  condenser  performance  is 
to  more  accurately  predict  condenser  characteristics  such  that  the  designer  can  develop 
reliable,  smaller,  and  hence  less  costly  condensers  which  meet  the  same  thermal  duty 
requirements.  Analyses  of  single  tubes  allow  isolation  of  individual  factors  which 
affect  condensation  without  the  complexities  which  result  from  vapor  flow  over  tube 
bundles  or  condensate  inundation.  Since  the  pioneering  work  of  Nusselt  [5],  a 
significant  amount  of  theoretical  work  has  been  done  analyzing  laminar  film 
condensation  on  horizontal  circular  tubes  (see  Rose  [6]),  the  results  of  which  have 
greatly  contributed  to  the  ability  to  reasonably  predict  single  tube  heat  transfer 
performance. 

Over  the  last  two  decades,  major  effort  has  been  expended  to  study  effects  which 
enhance  laminar  film  condensation  heat  transfer  over  that  obtained  from  a  plain 
circular  tube.  The  majority  of  these  techniques  have  focussed  on  controlling  the 
thickness  of  the  condensate  film  as  this  is  the  major  resistance  to  heat  transfer.  The 
condensate  film  thickness  generally  increases  with  streamwise  distance  from  the  top  of 
the  tube  and  is  dependent  on  both  the  rate  of  condensation  and  the  interfacial  shear 
between  the  vapor  boundary- layer  and  condensate  film.  Additionally,  at  high  vapor 
velocities,  vapor  boundary- layer  separation  occurs  which  results  in  a  rapid  thickening 
of  the  condensate  film  downstream  of  the  separation  point.  Enhancement  techniques 
include,  but  are  not  limited  to,  extended  surfaces  (fins),  profiled  tube  surfaces  (roped 
or  corrugated  tubes),  and  non- circular  tube  geometries.  Marto  [7],  Bergles  et  al.  [8]  and 
Webb  [9]  provide  comprehensive  reviews  on  such  enhancement  techniques.  These 
techniques  enhance  heat  transfer  through  an  increase  in  surface  area  to  volume  ratio 
and/or   the  use  of  surface   tension   to  thin  the  condensate    film. 


Non-circular  geometries  contribute  to  a  thinning  of  the  condensate  film  by  placing 
more  of  the  surface  in  line  with  the  direction  of  gravity;  thus,  making  the  average 
streamwise  gravity  component  in  the  momentum  balance  of  the  condensate  film  larger 
than  that  obtained  from  a  circular  geometry.  Dhir  and  Lienhard  [10]  applied  a  simple 
Nusselt  type  analysis  to  an  arbitrary  plane  or  axisymmetric  body  in  which  the 
streamwise  gravity  component  varied  with  streamwise  length.  The  results  of  this 
analysis  gave  expressions  for  condensate  film  thickness  and  Nusselt  number  which  were 
identical  to  those  derived  by  Nusselt  except  that  the  gravitational  constant,  g,  is 
replaced   by  an  effective    constant,   gc(i,  which  is  defined   by 

_     x  (g  R)4'3 

8eff"  ry^  (I"1} 

Jo 

where  x  is  the  streamwise  distance  from  the  leading  edge  of  the  surface  and  R  is  the 
radius  of  curvature  of  the  axisymmetric  body.  These  expressions  can  only  be  applied 
to  systems  for  which  curvilinear  coordinates  are  applicable,  i.e.,  for  bodies  in  which  the 
radii  of  curvature    are  much  greater   than  the  film  thickness. 

Shklover  and  co-workers  [11,12]  analyzed  a  horizontal   cylinder  with  a  logarithmic 
spiral  type  of  surface   curvature    (see  Figure  I- 1)  such  that  more  of  the  body  surface  is 
aligned     in    the    direction     of    gravity     and    the    radius    of    curvature     is    continuously 
increasing  along  the  streamwise   direction.   They  manufactured    this  tube  by  "mechanical 
deformation     of    a   circular     tube."      The    effect    of    this    profile    was    to   increase    the 
streamwise    effective    gravity    component  over  that  found   for  a  circular    cylinders    and 
to  utilize    surface    tension    effects    on  the  pressure    gradient    to  increase    the  mean   film 
velocity,  u(x,y),  in  the  streamwise   direction.    Their  analysis   considered    viscous,  gravity 
and    surface    tension    forces    in   the    momentum    balance    of   the    condensate     film    and 
neglected   vapor  shear  at  the  condensate/vapor    interface.     These  effects  resulted   in  a 
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Figure   1-1.  Logarithmic    Spiral   Tube  from   Reference    [12]. 


thinning  of  the  condensate  film  and  an  increase  of  20-30%  in  the  overall  heat-transfer 
coefficient  when  compared  to  a  circular  tube  of  the  same  surface  area.  These  results 
agreed  with  experimental  data.  In  analyzing  the  individual  effects  of  gravity  and 
surface  tension,  it  was  determined  that  surface  tension  was  most  significant  over  the 
top  portion  of  the  tube  where  dR/dx  is  largest  and  accounted  for  10-15%  of  this 
increase. 

An  elliptical  tube  with  major  axis  aligned  with  the  direction  of  gravity  has  a 
larger  effective  gravity  than  a  circular  cylinder,  as  well  as  an  increasing  radius  of 
curvature  over  the  top  half  of  the  tube.  In  addition,  they  are  more  practical  from  a 
manufacturing  standpoint  than  the  non-circular  tubes  analyzed  above.  Vapor 
boundary-layer  separation  should  also  occur  at  a  point  further  downstream  as  compared 
to  a  circular  tube  and  thus  the  rapid  thickening  of  the  condensate  film  would  e 
delayed.  Methods  of  manufacturing  elliptical  tubes,  and  associated  problems  will  be 
discussed    in  Chapters   IV  and  VI. 


Wallis  [13]  conducted   flow  visualizations     (Figure  I-2)of  water  around  circular   and 
elliptical    tube  bundles.     These  photographs    provide   a  qualitative     feel  for  the  flow 


Figure   1-2.   Flow   Visualization,     Water  over  Circular    and  Elliptical    Tube  Banks, 
from  Reference    [13]. 


characteristics    around   each  type  of  tube  bundle.    In  examining    the  photographs,    it  can 
be  concluded    that  the  separation    point  for  the  elliptical    tube  bank  shifts   downstream 
compared    to    the    circular     tube    bank    and    results    in    smaller     wake    regions    behind 
individual    tubes.    These  effects  should  contribute    to  better  heat  transfer    performance 
for  the  elliptical    tube  bank.    Experimental    studies  by  Joyner  and  Palmer  [14]  show  that 
single   phase   flow  resistance    and  pressure   drop  are  significantly     smaller   for  elliptical 
tube   bundles   as  compared    to  circular    tube   bundles.     In  power   condensers,    shell-side 
(vapor)  pressure  drop  has  an  effect  on  the  overall  efficiency    of  the  power  plant  and  on 
the  thermal   driving    force  for  condensation    (Tsat-Twan).  Minimizing   the  vapor  pressure 
drop  improves    power  plant   efficiency    and   maintains    a  relatively    constant    Tsat  (for  a 


large   pressure   drop,  Tsat  is  reduced   as  vapor   flows  through   the  condenser    resulting    in 
less  thermal   driving    force). 

Studies  involving  single-phase  heat  exchangers  and  condenser-evaporators  indicate 
that  elliptical  tube  heat  transfer  performance  is  superior  to  a  circular  tube  of 
comparable  surface  area.  Ota  and  Nishiyama  [15]  conducted  an  experimental  investi- 
gation of  single  phase  forced  convection  heat  transfer  of  air  over  an  elliptical  cylinder 
of  minor-to-major  axis  ratio  of  1:3  at  various  angles  of  attack.  They  concluded  that 
elliptical  cylinders  gave  improved  heat  transfer  performance  over  circular  cylinders  at 
all  angles  of  attack  (as  a  result  of  fluid  turbulence)  in  the  range  of  Reynolds  numbers 
studied  (8000  -  79000).  Moalem  and  Sideman  [16]  conducted  a  theoretical  analysis  of  a 
horizontal  elliptical  tube  in  a  condenser-evaporator  used  in  desalinization  plants.  Their 
study  showed  a  10-20%  enhancement  for  elliptical  tubes  as  compared  to  circular  tubes 
with  the  maximum  enhancement  achieved  at  a  minor-to-major  axis  ratio  (aspect  ratio) 
of  1:4  where  the  major  axis  is  aligned  with  the  vertical.  Huang  and  Mayinger  [17] 
conducted  an  experimental  free  convection  heat  transfer  study  around  elliptical  tubes 
and  found  optimum  improvement  in  heat  exchanger  performance  for  tubes  with  major 
axis  aligned  vertically  and  with  an  aspect  ratio  of  1:2.  Merker  and  Ba'hr  [18]  used  an 
analogy  between  momentum  and  heat  transport  to  derive  a  semi-empirical  relation  from 
which  to  determine  mass  transfer  rate  from  the  fluid  boundary- layer  into  the  free 
stream.  The  results  of  their  study  showed  improved  heat  transfer  performance  for  heat 
exchangers  in  which  the  elliptical  tubes  were  spaced  wider  in  the  longitudinal  direction 
and  more  compact  in  the  transverse  direction.  In  all  these  studies  it  appears  that  an 
elliptical    geometry   improves  the  heat  transfer    performance    of  a  heat  exchanger. 


B.        SIGNIFICANT      THEORETICAL      STUDIES     ON     HORIZONTAL      CIRCULAR 
TUBES 

The  purpose  of  this  survey  is  to  provide  a  historical  perspective  to  the 
understanding  of  laminar  film  condensation  on  a  horizontal  circular  cylinder.  Only 
those  works  which  have  made  significant  advances  in  this  understanding  are  considered. 
The  methods  of  some  of  these  analyses  will  later  be  applied  to  formulate  and  solve  the 
governing  equations  for  the  case  of  the  horizontal  elliptical  cylinder.  Figure  1-3 
provides   the  geometrical    layout  and  coordinate    systems  used  in  these  analyses. 


Figure     1-3.       Geometry     for     Film     Condensation      on    a 
Horizontal    Circular    Cylinder. 


Nusselt  [5]  used  a  simple  momentum  and  energy  balance  to  determine  the  heat 
transfer  properties  for  condensation  on  flat  plates  and  horizontal  circular  cylinders. 
In  simplifying    these  physical    laws  for  solution,   he  made  the  following    assumptions: 


(1)  The  tube  wall  temperature    is  constant. 

(2)  The  thermophysical    properties   of  the  fluid   are  constant. 

(3)  The  temperature    at  the  film/vapor    interface    is  Tsat. 

(4)  The  condensate    film   thickness    is  small  compared    to  the  radius   of  the  tube. 

(5)  Heat  transfer  in  the  condensate  film  is  one  dimensional  in  the  radial  direction 
which  implies  a  linear  temperature  distribution  when  the  film  is  very  thin. 
Convection    in  the  condensate    film   is  neglected. 

(6)  The  only  forces  acting  on  the  condensate  element  are  due  to  viscosity  and 
gravity.     Inertial    forces  are  neglected. 


(7)  Flow  in  the  condensate    film   is  laminar. 

(8)  The  vapor   is  quiescent    with  no  interfacial    shear. 

The  governing    equations    for  this  system  become 

1, — -  +  P[g  sin<t>  =  0 
dxL 


m  kfg    =    ^Tsat   "    T»aU> 


(1-2) 


where  m  is  the  condensate    mass  flux  rate  across  the  interface.     The  mean  heat  transfer 
is  given  by 

(I-3a) 


«Nu  -  0.943 


A,  g  ^  p,  sin<t> 
r\,  L  AT 


for  a  flat  plate  at  angle  <t>  to  the  horizontal    (not  valid   at  4>  =  0)  and 
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for  a  horizontal    circular    cylinder. 

Sparrow  and  Gregg  [19,20]  extended  the  Nusselt  [5]  analysis  by  applying  boundary 
layer  theory  to  the  condensate  film.  This  study  incorporated  convection  and  inertia  in 
the  energy  and  momentum  balance  while  continuing  to  neglect  interfacial  shear.  They 
determined  that  for  practical  engineering  fluids,  inertia  and  convection  had  negligible 
effect. 

Shekriladze  and  Gomelauri  [21]  included  interfacial  shear  by  considering 
momentum  transfer  across  the  interface  due  to  the  condensation  process.  They  used  an 
asymptotic  expression  for  the  interfacial  shear  based  on  an  infinite  condensation  rate 
given  by 

t«  ■  m(l/4  -  «,)  (1-4) 

where  U^  is  the  streamwise    velocity   of  the  vapor  at  the  outer  edge  of  the  boundary- 
layer.    It  was  also  assumed  that  U^»u6  such  that  the  film  velocity   may  be  neglected  and 


that  inertia  and  convection  effects  in  the  condensate  film  could  also  be  neglected.  As 
noted  by  Rose  [6],  the  validity  of  neglecting  ufi  requires  that  X|AT/ti,hf  be  small  (large 
Prandtl  number).  For  forced  convection  heat  transfer  on  a  horizontal  circular  tube, 
Shekriladze    and  Gomelauri    determined    the  following    expression: 

Nu  =  0V/&  (I"5) 

where 

&  -  ^^  d-6) 

is  the  two- phase   Reynolds    number.     When  the  effect   of  gravity    cannot   be  neglected 
(mixed   convection)    the  mean  heat-transfer    coefficient    is  given  by 

Nu  =  0.64^^1  +\/l  +  1-69  F  ^"7^ 

where 

F  m  Vklf  (1-8) 

A,  AT  Ul 

and  measures  the  relative  effects  of  vapor  velocity  and  gravity.  It  should  be  noted  that 
based  on  potential  flow  outside  the  vapor  boundary  layer,  the  interfacial  shear  will 
always  be  positive  and  therefore  this  method  does  not  predict  vapor  separation  as  would 
occur  in  reality.  Based  on  a  minimum  separation  angle  of  82°  for  single  phase  flow  over 
a  circular  cylinder  without  suction  and  no  heat  transfer  after  separation,  Shekriladze 
and  Gomelauri  conservatively  estimated  that  the  heat- transfer  coefficient  would  be 
reduced  by  about  35%.  Actual  separation  occurs  at  a  point  further  downstream  and 
results  in  a  heat- transfer  coefficient  between  the  result  for  separation  at  82°  and  no 
separation.  These  expressions  have  been  found  to  be  reasonably  accurate  for  cases  of 
high  condensation    rate  and  are  useful   for  their  simplicity. 

Fujii   et  al.   [22]   conducted    an   analysis    of   mixed    convection    condensation     on 
horizontal    cylinders   using  the  two-phase  boundary-layer    equations.    Their  formulation 


neglected   inertia,   convection   and  pressure  gradient  effects  in  the  condensate    film.   The 
governing   equations    for  the  model  are  given  by 


dU       dV      _ 
—  +  —  =0 

dx        dy 

vdU  +  vdU  =  v  dV^  +  r^fy 

dx  ay  *  dx         pv  dy2 


(I-9a) 


for  the  vapor  boundary-layer    and 
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for  the  condensate    film.    The  compatibility    equations    at  the  interface    are  given  by 


(I-9c) 
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An  approximate  integral  solution  of  the  momentum  equation  for  flow  over  a  cylinder, 
with  suction  by  Truckenbrodt  [23],  was  modified  to  more  closely  agree  with  the 
numerical  solutions  of  Terril  [24].  This  method  enabled  the  determination  of  the 
interfacial  shear  and  subsequent  solution  of  the  governing  equations  for  the  condensate 
film.  This  technique  also  predicted  the  point  of  vapor  separation  and  subsequently,  a 
more  accurate  heat-transfer  coefficient.  For  locations  downstream  of  the  vapor 
boundary -layer  separation  point,  the  interfacial  shear  was  assumed  to  be  negligible  and 
simple  Nusselt  theory  used  to  determine  the  heat  transfer.  Local  and  mean  Nuss^  it 
numbers  were  numerically  determined  for  boundary  conditions  of  uniform  wall 
temperature    and   uniform    wall   heat   flux.     Their   results    were  in  fair   agreement    with 
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with    experimental     data.       In    addition     to    the    dimensionless     parameter,     F,    a    new 
dimensionless    parameter,    G,  was  defined    as 

AIA, 


G 


p,  n 


/ T'/  (I- 10) 


1/M    Pvlv 

which  is  a  measure  of  the  condensation  rate.  As  G-°°,  the  results  of  Fujii  et  al.  [22] 
agreed  with  the  previous  results  of  Shekriladze  and  Gomelauri  [21].  Lee  and  Rose  [25] 
noted  an  error  in  the  formulation  of  Fujii  et  al.  [22]  in  that  the  interfacial  shear  terms 
in  the  condensate  film  momentum  equation  should  have  been  divided  by  2.  They  re- 
analyzed the  system  of  equations  using  both  the  modified  and  unmodified  method  of 
Truckenbrodt  [23]  and  found  little  difference  whether  one  used  the  modified  or 
unmodified    method  of  Truckenbrodt. 

To  date,  Gaddis  [26]  has  conducted  the  most  comprehensive   study  of  condensation 
on  horizontal   circular  tubes  using  the  two- phase  boundary- layer  equations.    His  analysis 
neglected  surface   tension  in  the  momentum  equation   as  well  as  viscous  dissipation    and 
pressure    in  the  energy   equation    of  the   condensate    film.     He  determined    that   inertia 
effects   were  negligible    for  most  media  (with   the  exception   of  liquid    metals)   and  that 
convection    effects   were  somewhat   significant    for  viscous   liquids    with  Pr  »■  1,  similar 
to  the  conclusions   of  the  simpler  analyses   of  Sparrow  and  Gregg  [19,  20].  He  identified 
three  regimes  of  behavior    for  flowing   vapors.    For  low  vapor   Reynolds   numbers,  Rev, 
gravity  effects  dominate  the  heat  transfer   and  Nusselt's  analysis   adequately    represents 
the  heat  transfer    behavior.     For  high  Rey,  condensation    was  shear  controlled    and  the 
Nusselt  number  is  proportional    to  the  square   root  of  Rev  as  found  by  Shekriladze    and 
Gomelauri    [21].    The  intermediate    region,  in  which  gravity   and  vapor  shear  are  both 
significant,    results   in  no  simple  relation.     Gaddis  also  analyzed    several   cases  of  flow 
separation    (vapor  boundary- layer   and/or  condensate    film  separation).     For  low  Rev  (  = 
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10  ),  vapor  boundary-layer  separation  occurred  with  no  flow  reversal  in  the  liquid  film. 
This  result  is  expected  as  the  pressure  gradient  effect  on  the  condensate  film,  which 
tends  to  slow  the  condensate  flow  over  the  rear  half  of  the  tube,  is  a  function  of  the 
vapor  potential  velocity.  At  moderate  Rev  (~  103)  vapor  boundary-layer  separation  id 
condensate  flow  reversal  near  the  interface  occurred  with  a  resulting  decrease  in  the 
mean  fluid  velocity  and  rapid  thickening  of  the  condensate  film.  At  high  Rev  (~  105), 
no  vapor  boundary -layer  separation  occurred  but  the  condensate  film  separated  at  the 
wall  causing  a  sharp  reduction  in  mean  condensate  velocity  with  a  rapid  thickening  of 
the  film  as  a  result  of  a  significant,  adverse  pressure  gradient.  Fluid  separation  at  their 
respective  boundaries  is  defined  as  the  condition  at  which  the  shear  is  less  than  or  equal 
to  zero. 

diMarzo  and  Casarella  [27]  conducted  a  similar  analysis  using  a  more  general 
solution  technique  and  arrived  at  the  same  results  as  Gaddis  [26].  In  their  analysis  of 
the  flow  separation  phenomenon,  they  provided  guidance  on  determining  heat  transfer 
performance  once  separation  had  occurred.  At  low  Rev,  gravity  effects  dominate  and 
no  flow  reversals  or  thickening  of  the  condensate  film  occurred.  The  use  of  a  Nusselt 
type  analysis  with  no  interfacial  shear  would  seem  prudent  in  the  region  beyond  vapor 
separation.  At  moderate  to  large  /?ev,  pressure  gradient  effects  dominate.  The  adverse 
pressure  gradient  over  the  back  half  of  the  tube  causes  flow  reversal  and  a  rapid 
thickening  of  the  condensate  film.  In  this  case,  it  would  be  prudent  to  neglect  heat 
transfer  completely  beyond  the  separation  point  of  the  vapor-boundary  layer  or 
condensate    film. 

Rose  [28]  studied  the  effects  of  pressure  gradient  in  the  condensate  film.  In  this 
case,  the  pressure  gradient  on  the  condensate  film  is  due  to  the  pressure  gradient  of  the 
vapor,  as  determined    by  potential    theory,   which   is  impressed   on  the  condensate    film. 
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He  simplified   the  analysis  by  using  the  Shekriladze-Gomelauri    model  [21]  which  ignores 

inertia    and  convection    in  the  condensate    film  and  uses  the  asymptotic    expression    for 

the  interfacial     shear.     As  a  result   of  this  analysis,    a  dimensionless    parameter,    P,  was 

defined   by 

p  u    Py  K  "i  (Ml) 

P,  t-i  &T 

which  represents  the  pressure  gradient  effect  (Gaddis  [26]  had  a  similar  combination 
of  dimensionless  parameters  which  were  equivalent  to  P).  When  P  =  0,  the  governing 
equations  reduce  to  those  of  Shekriladze  and  Gomelauri  [21].  He  concluded  that  the 
effect  of  including  pressure  gradient  was  to  improve  heat  transfer  over  the  forward 
half  of  the  tube  since  the  pressure  gradient  over  this  region  is  favorable  (tends  to 
increase  the  mean  film  velocity  which  reduces  the  condensate  film  thickness).  As  a 
result  of  the  formulation,  he  found  that  for  cases  where  P  >  F  /  8,  a  critical  angle  at 
some  point  on  the  rear  of  the  tube  was  reached  where  d6  /d<{>  -  °°.  In  this  case,  it  was 
not  possible  to  obtain  a  solution  over  the  entire  tube.  It  was  postulated  that  this  critical 
angle  might  indicate  some  instability  followed  by  some  degree  of  waviness  or 
turbulence.  It  was  noted  that  this  critical  point  was  reached  prior  to  the  point  at  which 
the  condensate  film  separated.  For  conditions  which  permitted  solution  over  the  entire 
tube  (i.e.  P  <  F/&),  any  increase  in  heat  transfer  achieved  over  the  forward  part  of  the 
tube  was  balanced  by  a  decrease  in  heat  transfer  as  a  result  of  the  adverse  pressure 
gradient  over  the  back  half  of  the  tube  such  that  there  was  little  change  in  the  mean 
heat-  transfer  coefficient.  Since  the  pressure  gradient  due  to  potential  flow  of  vapor 
just  outside  the  vapor  boundary- layer    is  given  by 


^  ■  -p  ufi 

dx  y  *  dx 


■  -P,^  <'-12> 
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which  is  symmetric  around  the  surface  of  the  tube  from  front  to  back,  the  favorable 
effect  over  the  front  half  of  the  tube  is  exactly  offset  by  the  adverse  effect  over  the 
back  half  of  the  tube.  Referring  to  the  studies  conducted  by  Gaddis  [26]  and  diMarzo- 
Casarella  [27],  in  the  region  of  moderate  Rev,  condensate  film  separation  from  the  tube 
wall  does  not  occur  but  the  pressure  gradient  sufficiently  retards  the  film  such  that 
flow  reversals  occur  and  the  film  still  rapidly  thickens.  In  the  present  study,  the  effects 
of  pressure  gradient  using  the  analysis  of  Rose  [28]  resulted  in  a  condensate  velocity 
distribution  which  approached  zero  as  the  critical  angle  was  reached.  This  phenomenon 
may  be  a  result  of  the  model  used  in  formulating  the  problem.  Since  Shekriladze  and 
Gomelauri  [21]  used  an  asymptotic  expression  for  the  interfacial  shear,  negative 
velocities  cannot  exist  in  the  velocity  profile.  The  positive  velocity  distribution  is 
retarded  by  the  adverse  pressure  gradient.  In  the  case  of  P  >  F  /  8,  the  pressure  gradient 
is  large  enough  to  cause  complete  stoppage  of  film  flow  at  <f>c  and  the  film  thickness 
increases   rapidly. 

Krupiczka  [29]  examined  the  effects  of  surface  tension  due  to  film  curvature  on 
condensation  on  circular  cylinders.  He  used  a  simple  Nusselt  type  model  but  included 
surface  tension  in  the  momentum  equation.  In  his  development  he  did  not  assume  the 
film  thickness  was  much  le^s  than  the  radius  of  the  cylinder  until  after  the  inclusion 
of  the  surface  tension  term  to  account  for  the  curvature  of  the  film.  The  resulting 
equation  was  a  second  order  ordinary  differential  equation  which  required  two  initial 
conditions.  The  first  condition  was  given  by  the  symmetry  of  the  problem.  However, 
the  initial  thickness  was  not  obvious  and  was  arbitrarily  chosen  to  be  that  obtained 
from  the  Nusselt  model.  He  concluded  that  the  effect  of  surface  tension  was  small  on 
the  forward  part  of  the  tube  and  increases  in  significance  over  the  back  part  of  the  tube 
due  to  the  rapidly    changing    film  thickness.    This  significance    was  dependent   on  the 
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magnitude    of  dimensionless    parameter,  A,  given  by 

11/4 


A  =  -We 
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3r\  ,Ja 


d-13) 


P/  g  r3  Pr 

which  is  dependent  on  the  surface  tension  of  the  fluid  through  a  "modified"  Weber 
number  (We  =  a/(pgr~))  and  the  radius  of  the  cylinder,  r.  The  conclusion  was  that 
surface  tension  can  be  significant  for  small  diameter  tubes  and  wires  where,  as  r 
becomes  smaller,  the  film  thickness  becomes  relatively  significant  and  hence  the  film 
curvature  cannot  be  neglected.  For  most  practical  cases,  r-»h  and  therefore  the  surface 
tension  effect  due  to  film  curvature  can  be  neglected.  The  effect  of  assuming  the 
initial  film  thickness  to  be  that  given  by  Nusselt  [5]  was  checked  and  verified  to  have 
minimal  effect  on  the  mean  heat-transfer    coefficient. 

C.        SIGNIFICANT     THEORETICAL     STUDIES    ON    HORIZONTAL     ELLIPTICAL 
TUBES 

Some  theoretical  work  has  been  done  on  horizontal  elliptical  cylinders  in  a 
quiescent  pure  vapor  using  a  Nusselt  type  model.  The  physical  orientation  and 
coordinate    system  is  shown  in  Figure  1-4. 

Cheng  and  Tao  [30]  approximated  the  surface  of  an  ellipse  by  several  circular  arcs. 
They  analyzed  condensation  on  these  arcs  using  the  same  assumptions  as  Nusselt.  The 
ellipse  was  aligned  such  that  the  major  axis,  a,  was  aligned  with  the  direction  of 
gravity.  From  their  numerical  results  they  determined  that  the  heat-transfer 
coefficient  decreased  with  increasing  eccentricity,  k  (defined  as  the  ratio  of  the  minor 
to  major  axis).  Values  of  k  were  varied  from  that  of  a  vertical  flat  plate  (k  =  0)  to  that 
for  a  horizontal  circular  cylinder  (k  =  1).  The  mean  Nusselt  number  was  determined  by 
a  surface  area  weighted  average  of  the  mean  Nusselt  number  for  each  circular  arc  and 
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Figure   1-4.  Geometry  for  Film  Condensation    on  a  Horizontal    Elliptical    Cylinder. 


is  given  by 


Nu  =  0.2876 


Pr 
Ja 


i,4E^^/>^») 


(1-14) 


where 


Z,.(<|>)  =  -(sin6)-4/3f9(sine)1/3^8  +  C, 
3  •'flu 


(1-15) 


Cj  is  an  integral  constant  which  links  each  arc  and  /  is  the  number  of  arcs  used  in  the 
approximation.  In  the  practical  range  of  eccentricity  (0.3  -  0.6),  the  mean  heat  transfer 
coefficient  was  increased  by  10  to  18%  over  that  of  a  circular  tube  with  the  same 
surface  area.  Ali  and  McDonald  [31]  conducted  a  similar  type  analysis  as  Cheng  and 
Tao  without  the  circular  arc  approximation  as  a  first  estimate  for  condensation  on 
inclined    circular    tubes. 

Wang  et  al.  [32]  used  the  Nusselt  assumptions  to  analyze  condensation  on 
horizontal  elliptical  tube  for  which  the  major  axis  is  oriented  at  an  angle,  a.  wi  . 
respect  to  the  vertical    axis.    They  obtained   an  expression    for  the  mean  heat- transfer 
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coefficient    given  by 

r    .      2     .3      11/* 

[4      n,  A7      ^J     2*Vo  Jo  H/ 

where 

13/2 


Z(P)  =  sin^pP 
Jo 


a/> 


a2sin2(P-a)  +  62sin2(p-a) 


sin^pjp      .  (1-17) 


The  results  of  this  study  showed  that  the  maximum    mean  heat- transfer    coefficient    is 
obtained    when  the  major  axis  is  aligned    with  the  vertical    (a  =  0°  ).    These  theoretical 
results    were   validated    with   experimental    data   for  an  elliptical    tube   with  semi-minor 
axis  of  8  mm  and  semi-major  axis  of  22  mm.   Sheng  and  Cha'o  [33]  noted  that  the  mean 
heat-transfer   coefficient    was  incorrectly   determined  since  the  heat-transfer   coefficient 
is  averaged    over   the  surface    area   and   the   radial    distance,    r,  is  not  constant    for  an 
ellipse.    In  the  course  of  this  current  study,  the  problem  formulation    of  Wang  et  al.  [32] 
was  rerun  with  the  correct  expression    for  the  mean  heat- transfer    coefficient    (as  later 
derived   in  Chapter   II).    The  corrected   mean  heat- transfer    coefficient    of  the  elliptical 
tube  analyzed   by  Wang  et  al.  is  11.3%  larger  then  a  circular   tube  with  the  same  surface 
area.     Additionally,     it  should    be  noted    that    the   heat   transfer    performance    of   this 
elliptical    cylinder    was  better  than  the  circular   cylinder    for  angular   orientations    up  to 
a   ~  50°.     The   corrected    theoretical     results    are   approximately     5%  lower    than    those 
determined   by  Wang  et  al.  and  more  closely  agree  with  their  experimental    data.    Figure 
(1-5)  shows  the  heat  transfer    enhancement    for  elliptical    tubes  of  varying    eccentricity 
and  orientation    angle. 

Sheng  and  Cha'o  [33]  studied  the  effects  of  surface  tension  and  variable  wall 
temperature  on  condensation  on  a  horizontal  elliptical  tube.  The  remaining  assumptions 
of  the  analysis    were   the   same   as  those   above.     For   the   wall    temperature,    a  cosine 
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Figure   (1-3).   Elliptical    Tube  Mean  Nu  Compared   to  Circular    Tube  Mean  Nu  for 
varying    Eccentricity    and  Orientation    Angle. 

distribution    of  the  form 

Ar=  (j^.f^d  -Aco&Q)  (1-18) 

was  used  (which  has  been  shown  (Memory  and  Rose  [34])  to  be  in  good  agreement  with 
experimental  data)  where  the  value  of  A  depends  on  the  ratio  of  the  outside  to  inside 
heat- transfer  coefficient.  They  determined  that  variable  wall  temperature  affected 
local  values  but  not  the  mean  values  of  the  heat-transfer  coefficient.  For  the  surface 
tension  effect,  it  was  assumed  that  the  film  thickness  was  much  smaller  than  thu  -adius 
of  curvature  of  the  elliptical  surface  and  thus  the  surface  tension  effect  was  due  'ely 
to  the  curvature  of  the  tube  wall.  This  is  in  contrast  to  the  study  of  Krupiczka  [  or 
a  circular  tube  where  the  surface  tension  was  due  to  the  curvature  of  the  condc  te 
film.    Over  the  back  half  of  the  tube  the  radius  of  curvature    decreases   with    trean.      .e 
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distance  resulting  in  a  retarding  effect  on  film  flow.  As  the  flow  is  slowed,  the  film 
thickness  increases  rapidly  and  at  some  critical  angle  becomes  infinite  in  magnitude. 
This  phenomenon  is  similar  to  that  described  by  Rose  [28]  for  the  pressure  gradient 
effect.  The  surface  tension  effect  was  determined  to  be  significant  for  k  <  0.6  .  In 
comparison  to  a  circular  tube,  surface  tension  causes  a  favorable  pressure  gradient  over 
the  front  half  of  the  tube  resulting  in  a  greater  film  velocity,  thinner  film  thickness  and 
improved  heat  transfer.  It  has  an  opposite  effect  on  the  back  half  of  the  tube  which 
tends  to  negates  this  improvement.  The  gravity  component  of  the  momentum  balance 
is  the  driving  force  behind  the  enhancement  in  elliptical  tubes.  Consideration  of 
surface  tension  results  in  a  slight  decrease  in  the  mean  Nusselt  number  as  compared  to 
the  situation    in  which  the  surface    tension   is  neglected. 

A  potential    advantage    of  an  elliptical    tube  compared    to  a  circular    tube   is  the 
difference     in    vapor    flow    characteristics     as   a  result    of   a  better    streamlined     shape. 
Panday     [35]    developed     an    explicit     numerical     method     for    two    dimensional      film 
condensation     and   applied    it  to  the   case   of  downward     flowing    vapor    over   elliptical 
cylinders.      Convection    and   inertia    were   included    in  the  condensate    film   as  well   as 
surface  tension  and  pressure  gradient   (as  a  result  of  potential   flow  of  vapor  outside  the 
vapor  boundary-layer).     The  interfacial    shear  was  approximated    using  the  asymptotic 
expression   for  infinite   condensation    rate.   Several  errors  were  found  in  the  expressions 
for  surface    tension   and   pressure   gradient    as  a  result   of  an  incorrect    analysis    of  the 
differential    streamwise    length,  dx.   The  first  error  is  the  result  of  assuming   that  dx  = 
rd<f>.  This  relationship    assumes  that  the  radial   distance    from  the  centroid   of  the  ellipse 
is  constant   over  the  interval    of  the  parametric    angle.    This  fact  is  not  true  (as  will  be 
shown  in  Chapter  II)  and  results  in  an  error  whose  magnitude    is  dependent  on  the  step 
size  used  in  the  numerical    procedure.    The  second  error  involves   the  expression    for  r. 
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If  the  parametric    angle  is  measured    from  the  vertical   axis  of  the  ellipse   (as  appears   to 

be  the  case  in  Panday's   expression    for  vapor  velocity   potential),    then  r  should  be  given 

by 

x  =  b  sin<j> 

y  -  a  cos4>  (1-18) 


r  =  yx2  +  y2  =  aycos2$  +  /fc2sin2<J> 

Panday    used  an  expression    for  radial   distance   given  by 

r  =  asjsm2^  +  *Wcj>      .  t1'19) 

The  magnitude  of  the  resulting  error  is  dependent  on  the  eccentricity  of  the  ellipse. 
Additionally,  there  appears  to  be  an  error  in  the  gravity  component  of  the  momentum 
equation.  Panday  uses  a  body  force  given  by  pgsin0,  which  is  true  only  for  circular 
geometries  but  is  untrue  for  elliptical  geometries.  His  results  will  therefore  not  account 
for  the  improved  performance  which  are  obtained  by  placing  more  of  the  surface  in 
line  with  the  direction  of  gravity.  Panday's  conclusions  were  that  the  overall  heat 
transfer  was  reduced  for  elliptical  tubes  at  low  velocities  (as  compared  to  a  similar 
circular  tube)  due  to  a  rapid  thickening  of  the  film  at  the  stagnation  point  and  the 
overall  heat  transfer  is  increased  for  high  velocities  due  to  increased  interfacial  shear. 
These  conclusions  are  opposite  to  what  one  may  expect.  At  low  velocities,  the  film 
thickness  near  the  stagnation  point  for  elliptical  cylinders  is  thinner  than  the  equivalent 
circular  tube  as  a  result  of  increased  effective  gravity.  This  effect  results  in  improved 
heat  transfer.  Pressure  gradient  effects  due  to  potential  flow  and  interfacial  shear 
should  be  negligible  since  the  velocity  is  small.  At  high  velocities,  the  streamlined 
geometry  of  elliptical  cylinders  results  in  lower  interfacial  shear  and  pressure  gradient 
effects   as  compared   to  circular    cylinders.     These  effects   should   result  in  a  thickening 
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of   the    condensate     film    as   compared    to   circular    cylinders    and    hence    reduced    heat 
transfer    for  the  same  effective    diameter    and  free  stream    vapor   velocity. 

D.        OBJECTIVES    OF  THIS  THEORETICAL    STUDY 

The  results  of  previous  studies  indicate  that  elliptical  tubes  can  be  used  to 
increase  condensation  heat  transfer  in  low  vapor  velocity  condensers  as  compared  to 
circular  tubes  with  the  same  surface  area.  The  present  investigation  examines  the 
effects  of  vapor  shear,  pressure  gradient  and  surface  tension  in  laminar  film 
condensation  on  a  single  horizontal  elliptical  tube  with  its  major  axis  aligned  with 
gravity  and  the  free  stream  velocity.  Interfacial  shear  is  estimated  using  both  the 
simple  assumptions  of  Shekriladze  and  Gomelauri  [21]  and  the  more  complex  technique 
of  Fujii  et  al.  [22].  The  latter  case  calculates  the  angle  at  which  the  vapor  boundary- 
layer  separates  such  that  the  effect  of  reduced  drag  on  the  mean  heat-transfer 
coefficient    can  be  evaluated. 
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II.    THEORETICAL    DEVELOPMENT 

Following  the  theoretical  developments  used  for  the  case  of  laminar  film 
condensation  on  a  horizontal  circular  tube,  theoretical  models  are  developed  for  laminar 
film  condensation  on  a  horizontal  elliptical  tube  of  eccentricity,  k.  This  development 
starts  with  a  Nusselt  [5]  type  model  in  a  quiescent  vapor,  and  then  adds  forced  and 
mixed  convection  as  in  the  models  of  Shedkriladze  and  Gomelauri  [21]  and  pressure 
gradient  as  considered  by  Rose  [28].  The  pressure  gradient  takes  into  account  the 
effects  of  potential  flow  outside  the  vapor  boundary- layer  as  well  as  surface  tension. 
Finally,  a  model  is  developed  which  analyzes  the  vapor  boundary- layer  and  boundary- 
layer  separation  following  Fujii  et  al.  [22].  Pressure  gradient  and  surface  tension  are 
not  considered  in  this  model  due  to  the  complexities  introduced  by  the  use  of  the  two- 
phase  boundary- layer  equation.  Where  possible  these  elliptical  models  are  checked 
against  existing  theories  for  the  "limiting"  eccentricities  of  a  circular  cylinder  (£=0), 
vertical    flat  plate  (k=l)  and  horizontal    flat  plate  (&-°°). 

A.        FACTORS    RELATED   TO  ELLIPSE   GEOMETRY 

Consider  an  elliptical  cylinder  whose  cross-section  is  oriented  such  that  the  major- 
axis  is  aligned  with  the  vertical  as  shown  in  Figure  II- 1.  The  eccentricity,  k,  of  the 
ellipse  is  defined   by 

k=  *    .  (H-1) 

a 

where  a  and  b  are  the  semi-major  and  semi-minor  lengths,  respectively.  Functions 
related  to  the  geometry  of  the  ellipse  are  initially  developed  in  a  cartesian  coordinate 
system  (xl,yl)  whose  origin  coincides  with  the  centroid  of  the  ellipse.  This  is  then 
transformed    into  a  cylindrical    coordinate    system   (r,8)  where   r  is  the  radial    distance 
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from  the  ellipse   centroid    to  a  point  on  the  ellipse   surface   and  9  is  the  angle   measured 


from  the  vertical.     The  transformation     equations    are 

Xj  =  r  sin6 
y,   =  r  cos6 


(II -2) 


The  elliptical    surface  can  also  be  defined   by  a  parametric    angle,  <j>,  measured   from 
the  upper  semi- major  axis,  such  that  the  yt  coordinate    on  a  circle  of  radius  a  translates 


Figure  II- 1.  Geometry   for  Film  Condensation    on  a  Horizontal    Elliptical 
Cylinder    with  iMajor  Axis  Parallel    to  Gravity. 


to  the  y1  coordinate   on  the  ellipse  and  the  x±  coordinate   on  a  circle  of  radius  b  translates 

to  the  Xj  coordinate    on  the  ellipse.    The  transformation     equations    are  given  by 

xx  =  b  sin(J> 
y,  =  a  cos<}> 

Combining   Equations    (1-2)  and  (1-3)  results  in  a  relationship    between  r,8  and  <j>: 


(II -3) 
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sind>  =  —  sin6  (a) 
b 

cos<|>  =  -  cos6  (b)  (II -4) 

a 

tan6  ,  x 

tan<(>  =  (c)  . 


It  is  assumed   that  the  condensate    thickness,    6,  is  much  smaller   than  the  radius  of 
the  elliptical    surface.      Therefore,    the  equations    for   the   condensate    film   and   vapor 
boundary- layer   are  developed  using  an  orthogonal,   curvilinear    coordinate    system,  (x,y), 
where  x  is  the  streamwise    distance    along   the  elliptical    surface    and  y  is  the  distance 
normal  to  the  surface. 

1.        Radial   Distance,  r 

In  cartesian    coordinates,    the  surface   of  an  ellipse  is  given  by 

Xl  +  Zl  =  i  01-5) 

b2       a2 

The  radial   distance,   r,  can  be  determined    by 


*t*>  ■  f^J  ■  (II"6) 

Using    the    parametric     transformation      equations     (Equation  (II- 3)),    an    equivalent 
expression    for  the  radial   distance    is  given  in  terms  of  <J>  by: 

K4>)  =  a\/  cos2<t>  +  k2  sin2<J)     .  (II_7) 
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2.         Radius  of  Curvature    of  an  Ellipse,  R 

The  curvature,   k(x1),  and  the  radius  of  curvature,    R(xT),  of  a  function,   f(x1), 
is  given  by 

lf"(*,)| 


*(*,) 


[1  +  (/r/(*1))2]3/2  (II -8) 

R(xJ  =  IMxJ 


where  the  prime  denotes  differentiation     with  respect  to  xv  Solving  for  y^  as  a  function 
of  x1  using  Equation    (II-5),  determining    the  radius  of  curvature    using  Equation    (II- 8) 
and  transforming    the  results  to  parametric    coordinates    gives  the  radius  of  curvature    as 
a  function   of  4>: 

If®)  =  -  [sin2*  +  JkW^f2    .  (II -9) 

k 

3.         Streamwise    Length,  x 

Since  the  radial  distance  to  a  point  on  the  elliptical  surface  is  not  constant, 
the  streamwise  length  is  not  proportional  to  0  as  it  would  be  for  a  circle.  Additionally, 
as  eccentricity  varies,  the  streamwise  length  for  a  given  6  is  not  the  same.  To  obtain  an 
expression  for  x  in  terms  of  parametric  angle  <\>,  consider  first  a  point  on  the  ellipse 
surface  as  shown  in  Figure  II-2.  Moving  a  small  distance,  dx,  results  in  incremental 
changes,  dd  and  dr.   The  resulting   relationship    between  x,r  and  0  is  given  by: 

(dx)2  =  (dr)2  +  (r  J6)2       .  (n-10) 

Expressions    are   now   needed    for   dr  and   dd  in  terms   of  d<f>-     For  dr,  this   is  simply 

achieved   by  differentiating     Equation   (II-7): 

dr  =  l«  (*2-D  *W    m  („_!!) 

2  v/cos2*  +  £2sin24> 
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Figure    II- 2.     Enlarged     View    of   Point    on 
Ellipse. 


For  d8,  taking   Equation    (II-4a),  substituting    Equation    (II- 7)  for  r  and  differentiating 
the  resultant    expression    results   in: 

cos6  dB  =  ^^ .  (H-12) 

[cos2<J>  S  khm^f2 

Combining    Equations    (II -4b)  and  (II- 7),  and  substituting    for  cos6  in  Equation    (11-12) 


results   in  an  expression    for  dd  in  terms  of  d(f>: 

k 


dB  = 


d$    . 


(11-13) 


cos24>  +  £2sirr<t> 

Substituting    Equations   (11-11)  and  (11-13)  into  Equation   (11-10)  gives  an  expression    for 
dx  as  a  function   of  d<f>: 


dx  -  a 


Jk2  f— (jfc2-  l)2sin22* 


\     cos2<J>  +  ifc2sin2<|> 


d$ 


(11-14) 


This  may  be  integrated    to  obtain    the  streamwise    length.     For  compactness,    a  function 
X(4>)  is  defined   by 


X(4>)  =  a 


k2  *  ±.(k*-V)W24 

4 

Nj      cos24>  +  A:2sin2<j> 


(H- 15) 
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In  later  developments,   a  characteristic    length  is  used  to  non-dimensionalize 

the  heat  transfer   parameters    in  the  model.    An  effective    diameter,   De,  is  defined   as  the 

diameter    of  a  circular    cylinder    having   the  same  surface    area  as  an  elliptical    cylinder 

and  will  be  used  as  a  characteristic     length.    The  effective    diameter    is  given  by 

*£>    =  2[\(<b)d<\>      .  (11-16) 

Jo 


Streamwise    length  is  non-dimensionalized     by: 

.        2x 
x     = 


7l£> 


(11-17) 


This    enables    direct    comparison    between    elliptical     and    circular    cylinders    since    it 
represents   an  equivalent    surface   area. 

4.         Streamwise    Gravity    Component,  gx 

The  component  of  gravity  in  the  streamwise  direction  is  tangent  to  the 
elliptical  surface  and  is  a  function  of  streamwise  location.  A  line  tangent  to  the 
surface,   defined   by  the  slope  of  the  surface,   is  obtained    from  Equation    (II-5), 

<*!?\  a        x\ 


dx. 


frT^ 


(11-18) 


A  unit  vector  in  the  streamwise    direction    is  then  given  by 

1 


T  = 


N 


,h2        2.       a\t 


hu2  2   f        a  f 

\jb    -  *i  i  -  -  Xj  j 


(11-19) 


The  streamwise  gravity  component  is  given  by  the  dot  product  of  the  gravity  vector  and 

the  tangent   vector, 

x, 


gx  =  g     T  =  g 


N 


bz 


(11-20) 
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Transforming    Equation    (11-20)  into  parametric    coordinates    results   in  the  expression: 

gx  =  g  fi(<\>) 

where     f^fy)  =  ± 

ysin24>  +  Jfc2cos2<|> 

5.  Potential    Flow  at  Ellipse  Surface,   U^ 

When  considering  forced  convection  due  to  the  flow  of  vapor  over  the 
condensate  film,  the  velocity  of  the  vapor  influences  the  condensate  film  thickness 
through  the  vapor  shear  at  the  film/vapor  interface.  The  vapor  velocity  at  the  interface 
will  be  dependent  on  the  vapor  velocity  outside  the  vapor  boundary- layer.  This  velocity 
is  determined  from  potential  flow  theory.  Assuming  that  the  film  thickness,  6,  and  the 
vapor  boundary- layer  thickness,  A,  are  much  smaller  than  the  radial  distance,  r,  of  the 
elliptical  surface  (and  therefore  may  be  neglected),  the  potential  flow,  LL,  about  an 
elliptical  surface  with  major  axis  aligned  with  the  vapor  free  stream  velocity,  U,,,,  is 
given  by 

<>♦  =  U-  ,     (l+k)  ■  ("-22) 

VI  +  it2cot?<fr 

Details   of  the  derivation    of  this  expression   are  provided   in  Appendix   (A). 

6.  Ellipse  with  Major  Axis  Perpendicular    to  Gravity 

For  the  case  where  the  ellipse  major  axis  is  perpendicular  to  the  vertical  as 
shown  in  Figure  II-3,the  previously  determined  equation  for  streamwise  component  of 
gravity  remains  the  same.  For  0  <  k  <  1,  the  major  axis  is  aligned  with  gravity.  For  1 
<  k  <  °°,  the  major  axis  is  perpendicular    to  gravity. 
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Figure     II-3.      Geometry     for    Condensation      on    Horizontal      Elliptical 
Cylinder    with  Major  Axis  Perpendicular    to  Gravity. 


B.        FREE    CONVECTION     CONDENSATION      ON    AN    ELLIPTICAL     CYLINDER 
(NUSSELT   [5]  TYPE  ANALYSIS) 

Consider  a  quiescent,  saturated  vapor  at  a  temperature  Tsat,  condensing  on  a 
horizontal  elliptical  cylinder  of  eccentricity,  k,  and  semi-major  axis,  a,  aligned  with  the 
direction  of  gravity  as  in  Figure  II-l.  The  same  simplifying  assumptions  are  made  as 
proposed  by  Nusselt  in  his  analysis  of  condensation  on  a  horizontal  circular  cylinder. 
These  assumptions    are: 

(1)  The  vapor  is  pure  and  quiescent. 

(2)  The  tube  wall  temperature    is  uniform   and  constant. 

(3)  The  thermophysical    properties   of  the  fluid  are  constant   and  evaluated    at  some 
given  reference    temperature. 

(4)  The  temperature    at  the  film/vapor    interface    is  Tsat. 

(5)  The  condensate    film  thickness    is  small  compared   to  the  radius  of  the  ellipse. 

(6)  Heat  transfer    in  the  condensate    film  is  one  dimensional    in  the  radial   direction 
providing    a  linear   temperature    distribution    across  the  film  (based  on  6-«R). 

(7)  The  only  forces  acting  on  the  condensate    film  are  viscous  and  gravity    forces. 

(8)  Flow   in  the   liquid    film    is  laminar    with    no  waviness    and   no  vapor   boundary 
layer  separation. 
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A  fluid   element   of  unit   width    is  analyzed    in  detail   and  shown   in  Figure   II-4. 


Figure      II-4.         Condensate       Film      Element      for      Free 
Convection. 


Conservation    of  mass  for  the  fluid   element  requires   that, 

d    f     r« 


Plf\  dy}  dx  -  Jp,/o*a  dy  +  -  jp^a  dy 

where  m  is  the  local  condensation    mass  flux  rate. 

Defining    the  mean   film   velocity    as 

um  =  —[  u  dy 
m       <5->o 


dx  +  m  dx  =  0 


Equation    (11-23)  simplifies    to 


m  -  -{pA 


(11-23) 


(11-24) 


(11-25) 


Conservation    of  momentum   for  the  film  element  is  a  balance  of  viscous  and  bodv 
forces   which   reduces   to 


dZU  r,x.s  a 

n,— -  +  p;  8  /i(<t>)  =  o 

dv1 


(11-26) 
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Here,  an  assumption    is  made  that  the  fluid   is  Newtonian,    i.e., 

t  =  n,  ^    -  (H-27) 

dy 

Conservation    of  energy   is  a  balance    of  latent    heat   from   condensation    and  heat 
conduction    through   the  condensate    film   which  reduces   to 


hfgm  dx  =  ^(7^-r^)  dx 


m  =  


(11-28) 


The  boundary    conditions    for  the  momentum   equation    are 

»,-o  -  ° 


*'  -  "•  ♦ 


=  0 

/y=8 


Integrating    the  momentum  equation  subject  to  the  boundary   conditions   results  in 
a  local  film  velocity   of 

"(y)  =  -^/1(4>)(8y-|vj  (H-30) 


1/ 
and  a  mean  film   velocity   of 


»„  =  — /i«»T    •  <n-31) 

1/  3 


Combining    the   continuity    equation,    energy   equation    and   expression    for   film    mean 
velocity   (Equations    (11-25), (11-28)  and  (11-31)),  results  in 

Using  the  relationship    between  dx  and  d<f>  from  Equations   (11-14)  and  (11-15)  and 


31 


carrying    out  the  differentiation     results   in 


hfg      i         3i,    X(d>)  {     "W'         /,WJ        dtj 


A  dimensionless    parameter,    Z},  is  defined    by 


1       5, 


(11-33) 


(11-34) 


where  B1  is  given  by 


*i  - 


p?  AA  * 


(11-35) 


Note  that  the  equivalent    parameter    used  by  Nusselt  [3]  multiplies    Equation    (11-35)  by 
a  factor  of  three  and  uses  radius  in  place  of  effective   diameter.    Substituting    Equations 
(11-34)  and  (11-35)  into  (11-33)  results  in  a  first  order  ordinary   differential    equation   for 
Z,: 


&\       4/ii(40  7         4  x(4>) 


*i  = 


(11-36) 


d<\>        3/W     l       D,/,») 


This  equation  can  be  solved  exactly  by  using  an  integrating  factor  which  would 
require  numerical  integration  or  can  be  solved  numerically  using  a  forward  stepping  or 
propagation  technique.  The  initial  condition  is  determined  based  on  the  symmetry  of 
the  problem  which  requires  that  (dZj/d<j>)A=0  =  0.  Solving  Equation  (11-36)  at  <f>  =  0 
results   in  the  initial    condition 

_    3  X(0)    _   3  a  k2 


J\  <J>=0 


Def[{0)  De 


(11-37) 
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The    non-dimensional      film     thickness,     6  ,    and     film     thickness,     6,    are    determined 

respectively    by 

6*(<t>)  =  Z,(<\>)m     and 

(11-38) 

6(<|>)  =  (  Z,(<|))  A,  )1/4    . 

The  local   heat-transfer    coefficient,     a,  and   mean   heat-transfer    coefficient    are  given 

respectively    by 

a((h)  =  and 

6(4>) 

r2*  ™  aa  (II"39) 


a  = 


/»' 


,  /"«(*)  x(4>)  d*   , 


flW. 


and  the  local  Nusselt  number  and  mean  Nusselt  number  by 

a(<|>)Z)  D 

Afa(<|>)  =  S  =  — !_      am3f 

X,  6(«|)) 

aZ)         2ml 

Afo  = =  - 1    — —  x(<l>)  <**      • 

Xt        Jo  5(<j>) 


(11-40) 


For  the  case  where  &  =  1  (circular    tube);    fj(<|>)=sin(t>,  f1'(<J>)=cos4>,  x(<t>)=a  and  De  = 
2a  and  Equation    (11-36)  reduces  to 

Si    +  i555*  z    =  —  (H-41) 

d<t>       3  sm4>     !       sine)) 

which  is  the  same  as  that  found  by  Nusselt  [5]  except  for  a  difference    in  the  definition 

of  Bv 

The  above  development    can  also  be  applied   to  a  horizontal    elliptical    tube  whose 

major  axis  is  perpendicular    to  the  direction    of  gravity   as  shown  in  Figure  II- 3.  In  this 

case,  b  is  greater   than  a,  i.e.  k  >  1. 
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C.        FORCED    CONVECTION    CONDENSATION     ON  AN  ELLIPTICAL    CYLINDER 
(SHEKRILADZE-GOMELAURI     [21]  TYPE  ANALYSIS) 

Consider    a   saturated     vapor    at    temperature     T    p    flowing     downward     over    a 
horizontal     elliptical     cylinder     with    free    stream     vapor    velocity,     Um.     The    elliptical 
cylinder    has    eccentricity,     k,  and    semi-major   axis,   a,  parallel     with    the    direction    of 
gravity    and  vapor   flow.    The  same  assumptions    as  for  the  case  of  free  convection    are 
used  here  with  the  exceptions: 

(1)    The  force  of  gravity    is  neglected. 

(3)  U,»uy=6. 

(4)  The    interfacial      shear    stress    is    approximated     by    an    asymptotic     expression 
assuming   an  infinite    condensation    rate, 

x6  =  m(U^u&)  =  mU^  (H-42) 


Conservation    of  mass  and  energy  are  the  same  as  previously  derived   (Equations 

(11-25)  and  (11-28)).    The  only   forces   acting   on  the  fluid    element  are   viscous    forces. 
Thus,  conservation    of  momentum   reduces  to 

Tl,  —t  '  0      ,  (""43) 


with  the  boundary    conditions, 

A  (11-44) 


du 


■  Ta  *  m  u* 


/>=6 

Integrating    the  momentum   equation    subject  to  the  boundary    conditions    results   in  a 
local  velocity   of 

u(y)    =     "L^*    y  (11-45) 
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and  a  mean  film  velocity 


um  =  ^*A     .  (H-46) 

Hi      2 


Using  the  expression  for  potential  flow  over  an  elliptical   cylinder  (Equation   (11-22))  and 
defining   a  velocity   function,    f2(4>),  such  that 

m)  =     l+k     ,  (n-47) 

the  mean  film  velocity   can  be  expressed   as 

«„  =  —/,<*)  f  •  (n-48) 

r\t  2 

Combining    the   continuity     equation,    energy    equation    and   expression    for   mean    film 
velocity   (Equations    (11-25), (11-28)  and  (11-48))  results  in 

X'AT  I  _  d  {pi  *i  ^  um         5)  (n_49) 

Using   the   relationship     between    dx   and   d<f>  from    Equations    (II- 14)  and   (II- 15),  and 
carrying   out  the  differentiation      results  in 

h*I  I  .  6  x-  Ar "•  _Lf,  ^  ,  #♦>  *|  .  (n-50) 

Defining   a  dimensionless    parameter,   Z->,  as 

Z2  =  *  (H-51) 

where 

s2  =  ^r    >  (II-52) 

p/  #. 


35 


Equation    (11-50)  becomes 


^1  +  2  ^  Z    =    4*(<t>)  (11-53) 


<ftj>         /2(<t>)  "2      £>e/2(4>) 

This  problem  is  a  first  order  ordinary  differential  equation  which  can  be  solved  using 
a  forward  stepping  numerical  integration  technique.  The  initial  condition  for  the 
problem  is  based  on  symmetry  which  requires  that  (dZ^/d^^Q  =  0.  Solving  Equation 
(H-53)  at  4>  =  0  results   in  the  initial    condition 

Z2„  =  ^W    .-**1    .  (,,-54) 

The  non-dimensional  film  thickness,  6  ,  and  film  thickness,  6,  are  determined 
respectively    by 


V  (II -55) 

6  =  yz2(<|>)  B2    . 

The  heat-transfer   coefficient    for  forced  convection   condensation    can  be  considered   in 


a  different    form.    Combining   Equations    (11-40), (11-51),  (11-52)  and  (11-55)  results  in 

1 

8*(4>) 


Nu  Re'm  =  — !—    ,     and 


Nu  Re'lr2  =  -P— —  x(<t>)<*t> 


(11-56) 

X(<t>)  <*t>      • 
'(♦) 

For  the  case  where  k  =  1  (circular    tube),  /2(4>)  =  2sin<|),  d*  =  ad$  and  Equation   (II- 

49)  reduces  to 

I  .  Llkft  sincD  «)    .  (H-57) 

6       ad(\>[    n, 

which  is  the  same  as  the  Shekriladze-Gomelauri     [21]  governing   equation    for  a  circular 
cylinder    with  no  body  forces. 
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D.        MIXED    CONVECTION     CONDENSATION     ON    AN    ELLIPTICAL     CYLINDER 
(SHEKRILADZE-GOMELAURI     [21]  TYPE  ANALYSIS) 

Consider  a  saturated  vapor  at  temperature  Tsat,  flowing  downward  over  a 
horizontal  elliptical  cylinder  with  free  stream  vapor  velocity,  U„.  The  elliptical 
cylinder  has  eccentricity,  A:,  and  semi-major  axis,  a,  which  is  parallel  with  the  direction 
of  gravity  and  vapor  flow.  The  forces  acting  on  a  film  element  are  viscous  and  body 
forces.  The  interfacial  shear  is  given  by  the  Shekriladze-Gomelauri  [21]  model  for 
infinite  condensation  rate  as  in  the  case  of  the  forced  convection  model  (Equation  (II- 
42)). 

Conservation     of  mass   and   energy   are   as  previously    derived.      Conservation    of 
momentum    for  the  condensate    film  reduces  to 

Vi—2   +  P,  S /,(<!>)  =0  (H-58) 

with  boundary    conditions, 

fJ.\  m    11  (H'59) 

Integrating    the  momentum   equation,   subject  to  the  boundary    conditions,    results 

in  the  local  film  velocity 

/ 

y 


u(y)  =  -^  /,(« 


5y- 


l   '       2  )  M/ 


Lmy  (H-60) 


and  mean  film  velocity 


/i(<t>)  —  +  /2(4>)  -    •  (11-61) 

1/  ■>  n;  2 


where   the  expression    for   potential    flow   over  an  ellipse    and   its  associated    velocity 
function   (Equations    (11-22)  and  (11-47))  have  been  used. 
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Combining  the  continuity   equation,  energy  equation  and  expression   for  mean  film 
velocity    (Equations    (11-25),  (11-28)  and  (II-61))  results   in 

h+z  i . -ihk  m  |  +  m^  m  a)  .        di-62) 

hfg       6        dx[i\,  3  T\thfg  2] 

Using  the  relationship    between  dx  and  d<f>,  and  algebraically    manipulating    Equation  (II- 
62)  results   in 

where  F  is  a  dimensionless    parameter    relating    free  and  forced   convection    and  6    is  a 
dimensionless    film  thickness   given  respectively    by 

1/  De  kfi  8 


F  = 

UlX,  A  T 

(11-64) 

o.  U  x    r^ 

5*  =  6 


p'  "-  =  A* 


Differentiating     Equation    (11-63)  results  in 

t:  ■  t^{2f«*>  6''  %  *  »A«  ^  ♦«»  £  *J*»  8'}  •    (II"65) 

5*        2  x(<t>)[  <*<{>  3  d<J>  J 

Regrouping    the  terms  shows  that   the  equation    is  a  first   order   ordinary    differential 
equation: 

{  2F/,(4»  6'3  +  m  6<}^  +  |  2F/;«t»  ^  +  )£<♦)}  =  1|$)     .  (H-66) 

The  initial  condition  for  this  problem  is  again  based  on  symmetry,  which  requires  that 
d6*/d<(>=0  at  <|>=0.  Inserting  this  condition  into  Equation  (11-66)  results  in  a  fourth  order 
expression    for  6 

|  F/fa)  6*4  +f2(0)  V1  -  2  ^)  =  o    ,  (11-67) 
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whose  root  of  interest   is  given  by 


6*(0)  = 

-^(0)  ♦ 

(^(0))2  +4{2/3F^(0)}{2   ^0)} 

2  (2/3  F/ftO)} 

in 


(11-68) 


3k 


1/2 


l+k\2+    16  Fa   _   Uk 
4F\\\  {   k  J  +     3  De  k 


Equation   (11-66)  can  be  solved  numerically    to  obtain  the  film  thickness.    The  other  heat 
transfer    parameters    are  determined    as  before  (Equations    (11-39)  and  (11-56)). 

For  the  case  where  F=0  and  £=1  (forced  convection  over  a  circular  tube),  Equation 
(11-63)  reduces  to  Equation  (11-57)  as  expected.  For  the  case  where  F  -  °°and  k  -\  (free 
convection    over  a  circular   tube),  Equation    (11-63)  becomes 

(11-69) 


1  D,      d   (2 

6*       2X(4>)   d*\3      /|W 


Upon   substitution     for   the   parameters    F  and   5     (Equation    (11-64)),  Equation    (11-69) 


becomes 


6        a  Xt  A7  i\f  d<\>{  3 


(11-70) 


which  is  the  same  as  Nusselt's  governing  equation   [5]  for  laminar   film  condensation    on 
a  horizontal    circular    cylinder. 

E.        MIXED    CONVECTION     CONDENSATION     ON    AN    ELLIPTICAL    CYLINDER 
WITH  SURFACE  TENSION   AND  PRESSURE   GRADIENT    EFFECTS 

Consider   a  pure  saturated    vapor  at  temperature    Tsat,  flowing   downward    over  a 
horizontal     elliptical     cylinder    with    free    stream    vapor    velocity,    U,,,.     The   elliptical 
cylinder    has   eccentricity,    k,  and   semi-major  axis,   a,  parallel    with    the   direction    of 
gravity   and  vapor  flow.    In  addition    to  the  assumptions    used  for  mixed  convection   in 


39 


the  last  section,  pressure  effects  in  the  condensate  film  are  considered  in  the 
development  of  the  governing  equations.  The  pressure  on  a  condensate  film  element 
will  be  a  combination  of  that  impressed  upon  the  element  by  the  potential  flow  of  the 
vapor  as  well  as  that  due  to  surface  tension.  A  fluid  element  is  shown  in  detail  in 
Figure   II-5. 


y 

~-4^ 

TU*U 

Figure     II-5.        Condensate      Film     Element      for     Mixed 
Convection    with  Pressure    Gradient. 


Conservation  of  mass  and  energy  are  as  previously  described.  Conservation  of 
momentum  for  the  film  element  is  a  balance  of  viscous,  pressure  and  body  forces  which 
reduces  to 

(H-71) 


with  boundary    conditions 


d2u  ,,,,       dp       n 


V-o  =  ° 
du\  m  u± 


<H-6  ^ 


(11-72) 
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As  previously    stated,   the  film   pressure   is  a  combination    of  the  pressure   due  to 

potential    flow,  p^,  and  the  pressure   due  to  surface   tension,  pa.  Thus  the  total  pressure 

gradient    is  given  by 

^  =  *♦   +  ^£     .  (11-73) 

dx         dx         dx 

The  pressure  gradient   due  to  potential    flow  is  given  by 


dp*  _  _      udU±  _    -pvU^   dU^ 
dx  Pv     *  dx  x(<t>)      <*t> 


(11-74) 


Substituting    the  expressions    for  potential    flow  and  x(4>)  (Equation    (11-22)  and  (11-15)), 
into  Equation    (11-74)  results   in 

(11-75) 


*♦    .         Pv   U- 


dx 


la 


m) 


where 


/,(♦)  = 


q+ky  k2  sin2<t> 
[  sin24>  +  k2cos2<b  ] 


\ 


cos24>  +  fc2sinz<j) 
k1  +  -(JfcM)2  sin2^ 


(11-76) 


The  pressure  due  to  surface   tension   is  given  by 

a 


P.  = 


m) 


It  is  assumed   that  5«eR.   Thus,  the  pressure  gradient   is  given  by 


dp0 


a    dR 


dR 


dx  gi  dx         /?2  x((j))   d$ 


(11-77) 


(11-78) 


Substituting    Equations    (II-9)  and  (11-15)  into  Equation    (11-78)  results  in 

-dx---2^m 


(11-79) 
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where 


m)  = 


k  (1  -k2)  sin2(|) 
sin2<j)  +  fc2cos24>] 


5/2 


\ 


cos  <(>  +  it  sin  4> 
k2  +  -  (*2-l)2  sin^ 


(11-80) 


Integrating    the  momentum  equation,  subject  to  the  boundary  conditions   results  in 
a  local  film  velocity   of 


u(y)  = 
and  a  mean  film  velocity, 


P/  8  „ .  s        Pv  u~   3       a 


l-  /,(♦> 

n,  2a  n, 


2    n^2 


L  ) 


♦^=/,(4»y  ("-81) 


Tl/  2a  n,  2    n;  a* 


^+^i^/2(4))8    .        (H-82) 
3  21,,    '2W' 


Combining  the  continuity   equation,  energy  equation  and  expression  for  mean  film 
velocity   (Equations    (11-25),  (11-28)  and  (11-82)),  results  in 


hfr      6       dx\ 


Tl/  2a  n,  2   ,,  a2 


p.  X.  A7  I/, 

— /2(<J))6 


3 


(11-83) 


Using  the  relationship    between  dx  and  </#,  and  algebraically    manipulating    Equation  (II- 
83)  gives 


JL  =   D<  d 

6*      x(40<*t> 


*/,<♦)  +  ^^/3«t»  +  | 


/D^2 


^  a  )  Bo 


/4(<t>) 


3 


(11-84) 


+  |/2«I>)SJ 


where  F  and  5    are  as  previously    defined   in  Equation    (11-64).  P  and  Bo  are 
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dimensionless    parameters    defined    by: 


Bo  = 


Py  hfg  1/ 
P,  h  Ar 

Pi  8  D] 


(11-85) 


P  is  related  to  the  pressure  gradient  and  Bo,  the  Bond  number,  relates  the  inertia  effects 
to  surface  tension  effects.  Differentiating  Equation  (11-84)  results  in  the  following  first 
order  ordinary    differential     equation: 


2  X(<t»   _ 


De  (Def     F 

2  F/,(<t>)  +  -5  p/3(<|>)  *  3  -4:    -f  /«(♦) 

a  \  a  J    Bo 


6'' 


♦#♦>«•}£♦ 


2FAd>)  +  —  F^(<tt 


(11-86) 


+  3 


'D* 


v  "  / 


Bo 


y 


Am  s*2 


In  the  previous  development,  film  separation  and  instability  were  not  applicable 
since  no  effects  were  considered  which  could  retard  the  film  flow  and  (d5  /d <(>)-»<» only 
at  (|>c=tc.  However,  in  the  present  analysis,  the  decreasing  potential  flow  and  increasing 
radius  of  curvature  of  the  elliptical  cylinder  causes  the  pressure  gradient  to  have  a 
retarding  effect  on  the  condensate  film  over  the  back  half  of  the  tube.  At  the  point  of 
condensate  film  separation  from  the  tube  wall,  the  film  shear  stress  at  the  wall  is 
negative  which  implies  flow  reversal  (du/dy  <,  0  at  y=0).  Differentiating  Equation  (II- 
81)  and  solving  for  the  film  separation  criteria  (du/dy  =  0  at  y  =  0)  results  in  the 
following    condition    for  separation    to  occur: 

(11-87) 


D.  D,\   F 

2F/1(<J>)  +  -£  J»/3(4>)  +  3   -£   -£-  /4(<t>) 

a  l  a  I  Bo 


b*    +  2/2(<J>)  6'  s  0 
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Condensate   film  instability    occurs  when  the  rate  of  change  of  film  thickness  is  infinite. 
From  Equation    (11-86),  this  condition    occurs  when 


2F/1(<J.)  +  -S  P/3(4>)  +  3 
a 


(D* 


;  a  )  Bo 


/4(<t>) 


8'    +f2(4>)  8*  =0    . 


(11-88) 


From  a  comparison  of  Equations  (11-87)  and  (11-88),  it  is  evident  that  the  condition  of 
film  instability  will  occur  first  (has  a  smaller  magnitude  than  the  film  separation 
criteria  for  all  conditions)  and  is  therefore  the  limiting  condition.  Rose  [28]  suggests 
that  this  instability  may  manifest  itself  in  wavy  or  turbulent  flow,  or  may  signify  the 
detachment  of  the  liquid  film  from  the  tube  surface.  For  the  case  where  1/Bo  =  0 
(surface  tension  neglected),  and  using  Equations  (11-21), (11-47)  and  (11-76)  for /1(<f>),/2(<()) 
and  /3(<J>)  in  Equation  (11-88),  it  can  be  seen  that  the  condition  for  which  solutions  can 
be  determined    over  the  entire  tube  surface   occurs  when 

P<-±(-L)2F  (H-89) 


For  a  circular  tube  (k=l),  Equation  (11-89)  reduces  to  P<F/S,  which  was  that  found  by 
Rose  [28].  If  the  opposite  of  Equation  (11-89)  is  true,  then  the  pressure  gradient  effect 
is  dominant  enough  to  significantly  retard  the  flow  over  the  back  half  of  the  tube  and 
result  in  a  rapid  thickening  of  the  condensate  film.  Note  that  the  condition  for 
condensate  film  instability  is  a  function  of  vapor  velocity  (F)  since  the  pressure 
gradient   due  to  potential    flow  is  velocity   dependent. 

For  the  case  where  P  =  0,  and  using  Equation   (11-88),  the  condition    which  allows 
solution   of  the  momentum   equation    over  the  entire  tube  is  now  given  by 


Bo       3 


/  _  \ 


KD., 


2     l4 


l-*2 


(11-90) 
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Note  here  that  the  condition  for  film  instability  is  only  a  function  of  geometry. 
Therefore,  for  small  values  of  eccentricity  (approaching  a  vertical  flat  plate),  the 
retarding  effect  of  surface  tension  on  the  condensate  film  over  the  back  half  of  the 
tube  is  significant  and  can  result  in  a  rapid  thickening  of  the  condensate  film  at  <f>c. 
The  initial  condition  for  this  problem  is  again  based  on  symmetry,  which  requires 
that  (d5  /d<J))^=0  =  0.  Combining  this  condition  with  Equation  (11-86)  results  in  a  fourth 
order  expression    for  6  : 


(D.\    F 


2  F//(0)  ♦  -^  Pfi(0)  OM    f  /fa) 
la  \  a  )  Bo 


*1  +/j0)  6*2  -  2*(0)   ,0     (H-91) 
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whose  root  of  interest    is  determined    by 


A  =  - 
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2Ffi(0)  *  -5  P&0)  +  3  p    I-  /[(0) 

a  \  a  )  Bo 


c  =    2x(Q) 
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(11-92) 


~\ 

5*(0) 

^B2-4  A  C 

-  B 

2>4 

Equation  (11-86)  can  be  solved  numerically    to  obtain  the  condensate   film  thickness.    The 
other  heat  transfer    parameters   are  determined    as  before  (Equations    (11-39)  and  (11-56)). 
It  can  be  shown  that  this  formulation    reduces  to  the  models  developed  for  circular 
tubes.   Consider  the  case  where  k=l  and  1/Bo=0  (circular   tube  with  no  surface   tension); 
Equation    (11-84)  reduces  to: 

—  -  —  I  [  2  F  sin<|)  +  8  P  sin2<t>  ]  —  +  2  6*  sin*  \      . 
6*       d<b[  3 
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which  is  identical    to  that  developed   by  Rose  [28]  where 

1 


5*  =  6 


t  IT  '^ 


l2flT,,J 


_L  6.  (H-94) 

v/2 


For   the  case   where   fc=l,  P=0  and    1/Bo=0   (circular    tube,   no  surface    tension    and   no 
pressure  gradient),   Equation    (11-84)  reduces  to  the  Shekriladze-Gomelauri     [21]  model. 

F.  MIXED  CONVECTION  CONDENSATION  ON  AN  ELLIPTICAL  CYLINDER 
WITH  VAPOR  BOUNDARY  LAYER  SEPARATION  (FUJII  ET  AL.  [22]  TYPE 
ANALYSIS) 

The  mixed   convection    model   developed    by  Shekriladze     and   Gomelauri    [21]  has 
been  used  because  of  its  simplicity   and  ease  of  solution.    However,  since  the  interfacial 
shear  is  approximated    by  an  asymptotic    expression   based  on  an  infinite    condensation 
rate    and    uses    potential     flow    outside    the    vapor    boundary-layer     (which    is   always 
positive),   vapor  boundary-layer    separation    does  not  occur.    As  described   in  Chapter  I, 
Fujii   et  al.  [22]   modified    the   interfacial     shear    stress    expression    by  simultaneously 
solving  the  boundary- layer  equations   for  the  condensate  and  vapor  ensuring  compatibil- 
ity   at   the    condensate/ vapor     interface.       This    technique    allows    for   a   more    precise 
description     of    the    condensation     problem    under    forced    convection     conditions     and 
includes   vapor  boundary- layer    separation. 

Consider  a  pure  saturated  vapor  at  temperature  Tsat,  flowing  downward  over  a 
horizontal  elliptical  cylinder  with  free  stream  velocity,  UM.  The  elliptical  cylinder  is 
oriented  such  that  its  major  axis  is  aligned  with  the  direction  of  gravity  and  vapor  flow. 
Due  to  the  complexities  of  solving  the  two-phase  boundary-layer  equation,  it  is  assumed 
that  the  only  forces  acting  on  the  condensate  film  element  are  gravity  and  viscous 
forces.  Surface  tension  and  pressure  gradient  are  neglected  which  is  reasonable  since 
they    were    found    to  have    negligible    influence    under    most   conditions    (discussed    in 
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Chapter  IV).  Other  assumptions  are  the  same  as  in  Section  II. E.  The  vapor  is  assumed 
to  be  in  thermal  equilibrium  such  that  the  energy  equation  for  the  vapor  need  not  be 
considered  in  the  analysis.  The  only  forces  acting  on  an  element  in  the  vapor  boundary- 
layer  (see  Figure  II - 6)  are  inertia,  viscosity  and  pressure  gradient,  which  is  impressed 
upon  the  element  by  potential    flow  outside   the  vapor  boundary- layer.     The  tangential 


Figure   II-6.    Condensate    Film   and   Vapor   Elements    for  Mixed  Convection    using 
Fujii  et  al.  [22]  Type  Model. 


velocity  at  the  vapor/condensate  interface  is  much  smaller  than  the  potential  velocity 
outside  the  vapor  boundary-layer  and  can  therefore  be  assumed  to  be  negligible.  Based 
on  the  above  assumptions,  conservation  of  mass  and  momentum  for  the  vapor  boundary- 
layer  are  given  by 


dU       dV 


dx  dy 
dl 
dx         pv    dy 


=  0        (continuity)  (a) 


„  dU      r,  dU      „    du*       nv  &U  „    x  ,.v 

U  —  +  V  —  =  (/.  — -  + (momentum)  (b) 

dx  dy  fir  n       A.2 


(11-95) 
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For  the  condensate    film,  conservation    of  mass,  momentum   and  energy  are  given  by 

+  g  /j((j))  =0       (momentum)  (a) 

P;   dy2 

(11-96) 

m  =  "HP/  "m  6)  =     '    , ~  ,  ~         (energy  and  mass)  (b)      . 

where   conservation    of  mass  and  energy   have   been  combined   as  in  previous   sections. 
Boundary    conditions    for  this  system  of  equations   are  given  by 

",=o  =  Vo  =  ° 

du\        __  0  <"-97> 


/y=6  V  vy  )    6 

(11-98) 


where  A  is  the  vapor  boundary-layer  thickness.  The  compatibility  conditions  at  the 
interface    require   that 

.  ,r  */  (*sat~*waU' 

and      -p/,  .  =  — — =-=-      , 

where  the  last  compatibility  relation  is  the  condensation  mass  flux  rate  (i.e.  for  a  given 
surface  area,  the  amount  of  vapor  condensed  must  equal  the  increase  in  the  mass  of  the 
film). 

The  technique  for  solving  the  vapor  boundary-layer  momentum  equation  involves 
an  approximate  integral  technique  developed  by  Truckenbrodt  [23].  The  initial  problem 
development  follows  the  method  of  Pohlhausen  [36].  The  continuity  and  momentum 
equations  (Equations  (II-95a)  and  (II-95b)),  are  integrated  over  the  thickness  of  the 
boundary  layer  from  y  =  6  to  y  =  6  +  A  and  then  combined  to  eliminate  the  normal 
component  of  vapor  velocity,   V,  resulting    in  a  momentum   integral   equation 


1     d  t    r      \        A,    dUM        V.  t. 

_L  f_(y2  a,)  +  — l-  — *  =  -*    +  L_  (11-99) 

Ul   **   *     2'       U,    dx         U,       pv  ul 
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where  Ajand  At  are  the  vapor  boundary   layer  displacement    and  momentum  thicknesses 
given  respectively    by 


wrKh 


(11-100) 


and      A,  =   j         — 


dy     , 


Vr  is  the  normal  component  of  vapor  velocity  at  the  vapor/condensate    interface    due  to 
the  condensation    process  itself  and  is  sometimes  called  the  vapor  boundary-layer  suction 
velocity.    Differentiating     the  A->  term  in  Equation   (11-99)  and  rearranging    the  equation 
results   in: 


d*2  +   *z  M* 
dx         U+     dx 


2   *  -i 

A2y 


u, 


(11-101) 
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With  some  further   algebraic    manipulation     this  yields 
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(11-102) 


To  solve   this   expression,    the   vapor    velocity    distribution     across    the   boundary-layer 
thickness    must  be  determined.     This  distribution    is  a  function    of  a  pressure   gradient 


parameter,   k,  and  a  suction  parameter,   k15  given  by  Tuckenbrodt    [23]  as 


K   = 


*2   Pv    <Kfr 

n..       dx 


and     k,  =  - 


V*  A2  Pv 


(11-103) 


Truckenbrodt    [23]  defined   a  shear  function    f(x,K1)  and  a  shape  function   H(k,Kj)  as: 

(II- 104) 


T*  A,  A, 

— — -        and       H(k,k.)  =  — 

1,  U,  '         A2 


and  a  dimensionless    function,    Z,  given  by 

Z  =  /fc. 


A2 


v*./ 


(11-105) 
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Equation    (11-102)  becomes 


dZ         u 
n    ££  =  2—  [/-  (2+H)  k  -  k.1 

e  dx  U>.[  *J 


(11-106) 


Based  on  existing   theoretical    data  (Schlicting    [37]  and  Torda  [38]),  Truckenbrodt    [23] 

proposed  a  simple   linear   approximation     for  Equation    (II- 106)  given  by 

dZ  =  _.  r0441(1_2K)  _  6k  i      .  (11-107) 

'dx        U.[  v  J 


Substituting    Equations   (II- 103),  (11-105)  and  (II- 106)  into  Equation   (II- 107)  and  defining 
dimensionless    velocities    by 


UM  =  —4        and       V.  = 


f/ 


(11-108) 


results   in 


dZ 
d$ 
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i  *♦>     0.441  f  1  *2VtM^Z  1  -  6  -^^*  Z 
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(11-109) 


The  dimensionless   suction  velocity,  Vfi,  is  determined  by  transforming    the  compatibility 
relation    involving    the  condensate    mass  flux  rate  in  Equation    (11-98)  and  is  given  by 

-VbjRFv=°-  (11-110) 

6* 

where  G  is  a  dimensionless   parameter  which  is  proportional   to  the  condensation    rate  and 

is  defined   by 


G    = 


if  K     N  Pv  n 


p/  n, 


(ii-iii) 


Combining    Equations    (11-109)  and    (11-110),  the   governing    equation     for    the    vapor- 
boundary    layer  becomes 


&=±m     0.441 
d*        ^     De    \ 


1  -2-^/Z 
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(11-112) 
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The  initial    condition    for  Equation    (11-112)  comes  from  symmetry   and  is  given  by 

— ]       =  0      .  (11-113) 


Combining   Equations    (11-112)  and  (11-113)  results   in 
D.    (dDA 


x(0)  I  <*t> 


(v/Z(0)  )2  +  0.882  G  —(0)  y/zm  -  0.441  =0 
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(11-114) 


This  is  a  quadratic    equation    from   which  Z(0)  is  determined    to  be, 
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(11-115) 


The  interfacial     shear   is  determined    from  Equation    (11-104).    Truckenbrodt    [23] 

provides   a  simple  approximation     for  the  shear   function,   /(k,Kj),  given  by 

(11-116) 


where 


/  (K,K,)    =    3.22    fa   (Ka+K) 


k.  =  0.0682  +  0.174  k, 


k  and  Kj  are  modified    by  combining   Equations    (11-103)  and  (11-105): 


De      dUt>  G   r- 
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(II- 117) 
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Dimensionless    interfacial    shear  is  defined   by: 


(11-119) 


PvUl 


An  expression   for  dimensionless   shear  can  then  be  obtained  by  combining  Equations   (II- 
104),  (11-105)  and  (11-116): 
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(11-120) 
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The  governing  equations  for  the  condensate  film  are  solved  in  the  usual  manner. 
The  momentum  equation  (Equation  (II-96a)),  is  integrated  to  solve  for  condensate  film 
velocity, 


U    =    f.(<\>) 


by  - 


2) 


—  y 


(11-121) 


The  mean  velocity   is  determined    from  Equation    (11-24): 

um  =  /,(<p)  —  +  6 

"      n,  /lW  3      2r,; 


(11-122) 


Substituting    Equation    (11-122)  into  the  energy/continuity     equation    (Equation    (II-96b)) 


and  multiplying    both  sides  of  the  equation    by 


Pi  u~  xi  (Tsa,-Twau> 


(11-123) 


results   in 


J.  =  -Sl  <L If m  *1  ♦  I  1  T6  a-: 
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(11-124) 


Carrying    out  the  differentiation     of  Equation    (11-124)  and  rearranging    results   in  the 
ordinary    differential    equation, 


X(4>)  _  F  dfx  x.<  _    1  J,    ^5  5.3 
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which,  due  to  the  symmetry   of  the  problem,  has  the  initial    condition, 

=  0      . 
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Combining   Equations    (11-125)  and  (11-126)  results   in 


X(0)  _  I 

D.         3 


)     5'4(0)-^N     6*3(0)=0  (H-127) 

4=o  4   G  {  d<S>  4=0 


from  which  6  (0)  can  be  determined    by  numerical    methods. 

For  the  case  of  k  =  1  (circular  tube),  Equation  (11-125)  reduces  to  the  analysis  of 
Fujii  et  al.  [22]  as  corrected  by  Lee  and  Rose  [25].  The  differential  equation  for  the 
vapor  boundary- layer  (Equation  (11-112))  and  for  the  condensate  film  (Equation  (II- 
125))  are  solved  simultaneously  using  the  compatibility  relations  of  Equations  (11-117), 
(11-118)  and  (II- 120).  At  the  point  of  vapor  boundary- layer  separation,  the  interfacial 
shear  stress  becomes  negative.  Downstream  of  this  point,  it  is  assumed  that  the  shear 
stress  at  the  vapor/condensate  interface  is  negligible  and  a  simple  Nusselt  type  analysis 
is  used  for  the  remainder  of  the  elliptical  surface.  The  heat  transfer  parameters  are 
determined   as  before  (Equations    (11-39)  and  (11-56)). 
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III.    NUMERICAL   METHOD 

In  the  theoretical  development,  the  governing  equations  for  each  model  were 
either  reduced  to  a  single  first  order  ordinary  differential  equation  (ODE)  or  to  a 
simultaneous    solution   of  a  system  of  first  order  ODEs  of  the  form 

&=Mw)    .  (ni-i) 

ax 

Initial   conditions    were  determined    from  the  symmetry  of  the  problem  which  required 
that  both  the  slopes  of  the  condensate    film  and  vapor  boundary-layer    were  zero  at  <J>  = 
0.   These  systems  of  equations   are  solved  using  one  of  the  forward   stepping  numerical 
methods. 

A.        SOLUTION     TECHNIQUE     FOR    THE    NUSSELT     [5]    AND    SHEKRILADZE- 
GOMELAURI    [21]  ANALYSIS    METHODS 

Each  analysis  method  resulted  in  a  single  first  order  ODE  of  the  form  given  by 
Equation  (III - 1).  The  solution  of  this  equation  provides  the  condensate  film  thickness 
along  the  elliptical  tube  surface.  Many  numerical  methods  are  available  to  solve  this 
type  of  problem,  each  with  its  own  advantages  and  disadvantages.  For  this  study,  a 
fifth  order  Adams  method  predictor-corrector  algorithm  was  developed  from  Crandall 
[39].  The  Adams  methods  are  multi-step  techniques  which  use  finite  difference  type 
operators  incorporating  previously  determined  points  (hence  the  term  multi-step)  to 
determine  the  value  corresponding  to  the  next  step.  These  methods  provide  the  same 
accuracy  and  are  more  efficient  than  the  Runge-Kutta  (single  step)  methods  but  are 
restricted  to  a  fixed  step  size.  The  Adams  method  employed  in  this  study  uses  an 
Adams-Bashforth  explicit  method  [39]  to  predict  the  value  of  the  dependent  variable  at 
the  next  step  and  then  an  Adams- Moulton  implicit  method  [39]  to  correct  this  predicted 
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value.  This  method  is  not  self  starting  and  requires  determination  of  the  first  four 
points.  To  start-up  the  algorithm,  lower  order  Adams  methods  were  used  to  obtain  these 
four  points.    The  recurrence    formulae    for  this  algorithm    are  given  by 

ys  =y*., +  h  Ep*/"(w*)  (in -2) 

where  x  is  the  independent  variable,  y  is  the  dependent  variable,  h  is  the  step  size,  s  is 
the  step  number  for  which  y  is  to  be  determined,  k  are  the  points  determined  in 
previous   steps  and   Pk  are  the  coefficients    provided    in  Table  III- land  III-  2. 

TABLE  III- 1.  RECURRENCE  COEFFICIENTS  FOR  PREDICTOR  FORMULA 


FORMULA 

s-5              s-4 

Pk 
s-3 

s-2 

s-1 

TRUNCA- 
TION 
ERROR 

(1) 

1 

0(h2) 

(2) 

-1/2 

3/2 

0(h3) 

(3) 

5/12 

-16/12 

23/12 

0(h4) 

(4) 

-9/24 

37/24 

-59/24 

55/24 

0(h5) 

TABLE  III-2.  RECURRENCE  COEFFICIENTS  FOR  CORRECTOR  FORMULA 


FORMULA 

s-4              s-3 

Pk 

s-2 

s-1 

s 

TRUNCA- 
TION 
ERROR 

(1) 

1/2 

1/2 

0(h3) 

(2) 

-1/12 

8/12 

5/12 

0(h4) 

(3) 

1/24 

-5/24 

19/24 

9/24 

0(h5) 

(4) 

-19/720       106/720 

-264/720 

646/720 

251/720 

0(h6) 

The  algorithm    used  in  the  solution    for  mixed   convection    condensation    with  surface 
tension  and  pressure  gradient   (Section  II.E)  is  provided   in  Appendix   (B). 
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For  many  applications,  this  type  of  start-up  routine  would  normally  not  provide 
sufficient  accuracy.  However,  no  difficulties  were  encountered  using  this  numerical 
procedure  for  solving  the  asymptotic  interfacial  shear  stress  models.  Solution 
convergence  was  obtained  for  an  angle  step  size  of  0.1°.  The  difference  between  an 
angle  step  size  of  1.0°  and  0.1°  was  less  than   1%. 

B.        SOLUTION    TECHNIQUE    FOR  THE  FUJII    ET  AL.  [22]  ANALYSIS    METHOD 

Reduction  of  the  governing  equations  resulted  in  a  system  of  two  simultaneous 
ODEs,  one  for  the  momentum  thickness  of  the  vapor  boundary- layer  and  one  for  the 
condensate  film  thickness.  The  Adams  method  previously  described  did  not  provide 
sufficient  accuracy  to  solve  this  two- phase  boundary  layer  problem.  In  particular,  the 
system  of  equations  exhibited  a  stiffness  problem,  i.e.,  the  numerical  method  could  not 
accurately  solve  the  vapor  boundary-layer  ODE  over  the  first  several  steps  even  with 
an  angle  step  size  of  0.001°.  To  obtain  sufficient  accuracy,  a  much  smaller  step  size 
would  be  required  which  would  significantly  reduce  the  efficiency  of  the  computer 
algorithm.  Lee  [40]  indicated  that  similar  difficulties  had  been  encountered  when  Lee 
and  Rose  [25]  solved  the  same  equations  for  a  horizontal  circular  cylinder.  Lee  and 
Rose  [25]  used  a  Runge-Kutta  method  which  permitted  variable  step  size.  Thus,  the  step 
size  could  be  reduced  sufficiently  to  obtain  the  required  accuracy  for  the  vapor 
boundary-layer  ODE  over  the  front  of  the  tube  and  then  increased  over  the  remainder 
of  the  tube  to  improve  the  efficiency  of  the  algorithm.  Utilizing  this  numerical, 
procedure,  they  were  able  to  solve  the  system  of  equations  over  most  of  the  range  of 
dimensionless  parameters  (F  and  G).  However,  in  some  cases  (specifically  large  values 
of  F  and  G),  the  system  of  equations    were  still  too  stiff   for  solution. 

To  overcome  the  problem  of  stiffness,  a  numerical  problem  solver  (IVPAG)  from 
the  International     Mathematical    and  Statistical     Library    (IMSL)  [41]  was   used.     This 
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algorithm    uses  Gear's  stiff   method  [42]  to  efficiently    solve  systems  in  which  stiffness 
is  a  problem.    However,  as  in  the  study  of  Lee  and  Rose  [25],  this  method  could  still  not 
solve  the  two-phase  boundary-layer    equations    over  the  entire   range  of  dimensionless 
parameters    considered,    but  did  cover  the  practical    ranges  of  the  parameters. 

It  was  originally  conceived  that  the  computer  algorithm  would  progress  around 
the  tube  surface  until  the  criteria  for  vapor  boundary- layer  separation  was  reached 
(dx6/d4>  <  0).  At  this  point,  the  model  would  shift  from  solution  of  the  two-phase 
boundary-layer  equations  to  the  simple  Nusselt  model.  However,  the  numerical  solver 
encountered  difficulties  in  solving  the  system  of  equations  as  it  approached  the  point 
of  vapor  boundary-layer  separation.  Therefore,  the  solution  method  was  changed  to  a 
two  step  procedure.  In  the  first  step,  the  system  of  equations  were  solved  until  the 
computer  algorithm  indicated  a  problem  had  been  reached  (usually  the  algorithm 
indicated  a  stiffness  problem  had  been  encountered).  The  interfacial  shear  stress  data 
was  analyzed  to  verify  that  this  problem  was  the  result  of  rapidly  changing  conditions 
associated  with  vapor  boundary- layer  separation.  Therefore,  to  facilitate  the  solution, 
vapor  boundary- layer  separation  was  assumed  to  have  occurred  if  the  interfacial  shear 
stress  at  the  last  point  obtainable  was  less  than  ten  percent  of  the  maximum  interfacial 
shear  stress  observed  over  the  surface.  At  this  point,  the  problem  solver  was  within  one 
to  two  step  sizes  of  the  actual  separation  point  (for  a  step  size  of  0.1°,  the  computed 
separation  point  was  therefore  within  0.2°  of  the  actual  separation  point).  The  second 
step  involved  re-running  the  analysis  using  the  computed  separation  point  calculated 
in  the  first  step  as  the  point  at  which  the  two- phase  model  was  switched  to  the  simple 
Nusselt  model.  A  converged  solution  was  obtained  for  an  angle  step  size  of  0.1°.  The 
algorithm    for  solving  the  two-phase  model  is  provided   in  Appendix  (C). 


57 


IV.    RESULTS   AND  DISCUSSION 

When   comparing    solutions    for   an  elliptical    tube    with    those   of  a  circular    tube, 
equivalent    surface    areas    have    been   used.     The  effective    diameter    of  the  equivalent 
circular   cylinder  is  given  by  Equation   (11-16).  The  dimensionless   streamwise   length  x  , 
defined  as  the  ratio  of  the  streamwise    length,  x,  to  the  half  perimeter  length  and  given 
by 

x*  =  A£_     ,  (iv-i) 

71  De 

is  used  when  comparing    local  values  of  film  thickness   and  heat- transfer    coefficients. 
Its  use  enables  direct  comparison  between  elliptical   and  circular  tubes  since  it  represents 
an    equivalent     fraction     of    the    total    perimeter     from    the    top    of    the    tube    to    the 
dimensionless    point. 

A  practical  range  of  eccentricities  for  an  elliptical  tube  with  major  axis  aligned 
vertically  is  0.3  <  k  <  0.6  .  The  lower  limit  of  0.3  is  based  on  discussions  with  a  tube 
manufacturer  (Reference  [43])  in  which  the  elliptical  tube  is  formed  by  pressing  a 
circular  tube  with  a  roller  assembly  to  obtain  the  proper  major  and  minor  axis 
dimensions.  Roller  contact  points  would  be  evenly  spaced  to  maintain  the  elliptical 
curvature  and  to  prevent  an  overly  flat  surface.  The  upper  limit  is  chosen  as  the 
eccentricity  in  which  noticeable  effects  were  observed  in  the  heat  transfer 
characteristics    (as  discussed   in  this  chapter). 

A.        EFFECT  OF  GRAVITY 

As  previously  noted  (Section  I. A),  some  improvement  in  condensation  heat  transfer 
is  expected  for  an  elliptical  tube  in  quiescent  vapor  compared  to  a  circular  one  since 
more  of  the  surface  is  aligned  with  gravity.    By  increasing   the  "effect"  of  gravity  on  the 
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condensate  film,  the  mean  condensate  velocity  at  a  given  streamwise  location  is 
increased  which  results  in  a  thinning  of  the  condensate  film  and  increased  heat 
transfer.  The  gravity  function,  f^cf)),  defined  in  Equation  (11-21)  is  shown  in  Figure  IV- 
1. 
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Figure  IV- 1.   Gravity    Function,    f^^)  vs.  Streamwise    Length,  x 


As  the  elliptical  tube  approaches  a  flat  plate  configuration  (k  -  0),  the  effect  of  gravity 
increases  to  sweep  the  condensate  film  along  the  tube  surface.  The  effect  of  gravity  on 
the  local  film  thickness  is  shown  in  Figure  IV-2for  two  eccentricities.  As  can  be  seen, 
the  condensate  film  is  thinned  over  the  front  and  rear  portion  of  the  tube  compared  to 
a  circular  tube  but  is  slightly  thicker  in  the  middle  region.  The  thickness  of  the  film 
is  controlled  by  the  rate  of  condensation  (which  is  dependent  on  condensate  film 
thickness)  and  the  condensate  velocity  (which  is  also  dependent  on  condensate  film 
thickness  and  gravity).    Near  the  top  of  the  elliptical   tube,  the  larger  gravity  component 
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increases    the   condensate     velocity    relative    to   that    of   a   circular    tube    resulting    in  a 
thinner     condensate      film.       This    thinner     film,     however,     results     in    an    increased 
condensation     rate    which    tends    to   thicken    the    film    further    downstream.      Since    the 
elliptical    tube  has  a  greater  condensation    rate  near  the  top  of  the  tube,  the  film  thickens 
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Figure   IV-2.   Local  Film  Thickness    for  Varying    Eccentricity    (k). 

more  rapidly  compared  to  that  on  a  circular  tube  causing  the  relatively  thicker  film  in 
the  middle  region  as  seen  in  the  figure.  At  90°,  the  gravity  component  is  the  same  as  on 
a  circular  tube  and  so  has  no  relative  thinning  effect.  Over  the  rear  half  of  the  tube, 
the  reduced  condensation  rate  (as  a  result  of  this  thicker  condensate  film)  and  increased 
velocity  due  to  the  larger  gravity  component  relative  to  a  circular  tube  results 
continued  thickening  of  the  film  on  the  elliptical  tube,  but  at  a  slower  rate  than  that 
on  a  circular  tube.  Eventually,  the  elliptical  tube  film  thickness  is  again  thinner  than 
on  a  circular    tube  for  the  same  dimensionless    streamwise    length.    The  condensate    film 
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thickness  is  therefore  determined  by  a  balance  of  the  effects  of  two  effects:  increased 
condensate  velocity  (which  tends  to  thin  the  film)  and  increased  condensation  rate 
(which   tends  to  thicken    the  film). 

The  overall  effect  on  the  mean  heat-transfer  coefficient  (given  as  a  mean  Nusselt 
number)  is  shown  in  Figure  IV-3.  To  compare  the  elliptical  tube  with  a  circular  tube, 
the  leading   coefficient    of  the  Nusselt  solution   (Equation    I-3b)  is  plotted  against 
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Figure  IV-3.   Mean  Nu  Leading   Coefficient    for  Varying   k. 

eccentricity.  Solutions  for  a  circular  cylinder  (k  =  1.0)  and  a  vertical  flat  plate  {k  -  0, 
L  =  2a)  agree  well  with  the  results  of  Nusselt  [5].  For  the  horizontal  flat  plate  (k  -  «), 
the  heat-transfer  coefficient  approaches  zero  as  would  be  predicted  by  the  Nusselt 
model.  For  the  practical  range  of  eccentricities  (0.3  <  k  <  0.6),  it  can  be  seen  that  the 
effect  of  placing  more  of  the  tube  surface  in  the  direction  of  gravity  is  to  increase  the 
mean  heat-transfer    coefficient    by  1%  for  k  -  0.6  and   13%  for  k  =  0.3.   Conversely,    by 
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placing  more  of  the  tube  surface  perpendicular  to  the  direction  of  gravity  (i.e., 
approaching  a  horizontal  flat  plate),  the  mean  heat-transfer  coefficient  decreases  as 
expected.     These  results   are  in  agreement    with  earlier   elliptical    tube  studies   [32,33]. 

B.         EFFECT  OF  VAPOR   VELOCITY 

1.         Asymptotic    Interfacial    Shear  Stress  Approximation 

The  streamlined  shape  of  the  elliptical  tube  has  an  effect  on  the  vapor  flow 
over  the  tube.  The  vapor  velocity  function,  f •,(<}>),  defined  by  Equation  (11-47)  is  shown 
in  Figure  IV-4  versus  dimensionless  streamwise  length  x  .  As  seen  in  the  figure,  the 
elliptical  tube  experiences  higher  vapor  velocities  than  a  circular  tube  at  the  front  and 
rear  of  the  tube  but  a  lower  vapor  velocity  in  the  central  region.    Using  a  balance  of  the 


Figure   IV-4.   Vapor  Velocity    Function,    f2(4>),  for  Varying   k. 
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factors  which  thicken  and  thin  the  condensate  film,  for  the  case  of  pure  forced 
convection  (F  -  0,  no  gravity  effect),  the  higher  relative  vapor  velocity  at  the  top  of  the 
elliptical  tube  results  in  a  larger  interfacial  shear  stress  and  hence  a  thinner  condensate 
film  and  larger  condensation  rate  compared  to  a  circular  tube.  With  more  condensate 
flowing  into  the  middle  section  of  the  elliptical  tube  and  a  lower  relative  vapor 
velocity,  the  condensate  film  thickens  more  rapidly  than  in  the  case  of  a  circular  tube. 
As  the  condensate  flows  over  the  rear  of  the  tube,  the  decreased  condensation  rate  (due 
to  the  thicker  film  in  the  middle  region)  and  increased  vapor  shear  results  in  thinner 
film  compared   to  the  circular    tube.    Figure  IV-5  shows  these  effects   in  detail   for  a 
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Figure   IV-5.    Effect   of  Vapor   Velocity    on  Local   Film  Thickness    for  a  Circular 
Tube  and  Elliptical    Tube  (k  =  0.6). 


circular    tube  (k  =  1)  and  an  elliptical    tube  (k  =  0.6)  using   the  asymptotic    interfacial 
shear  stress  (Shekriladze-Gomelauri)     approximation    and  the  two- phase  boundary   layer 
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shear  stress  (Fujii  el  al.)  approximation  (to  be  discussed  later).  Whereas  the  gravity 
effect  (f](4>))  described  in  Section  IV. A  was  always  larger  for  the  elliptical  tube  (see 
Figure  IV- 1),  the  vapor  shear  effect  f-,(<j>)  is  larger  or  smaller  depending  on  the  section 
of  the  tube  to  be  analvzed  (see  Figure  IV-4).  The  overall  effect  of  vapor  shear  usint  ne 
Shekriladze- Gomelauri  [21]  model  is  to  reduce  (by  about  2%  in  the  practical  range  of 
eccentricities)  the  elliptical  tube  mean  heat-transfer  coefficient  relative  to  an 
equivalent    circular    tube. 

In  the  mixed  convection  region,  the  mean  heat-transfer  coefficient  of  an  elliptical 
tube  is  increased  or  decreased  compared  to  a  circular  tube  depending  on  the  relative 
magnitudes  of  the  vapor  shear  and  gravity  effects.  This  relative  magnitude  is  measured 
by  dimensionless  parameter,  F.  For  large  F,  gravity  is  dominant  and  the  heat-transfer 
coefficient  trend  is  described  in  Section  IV. A  .  For  small  F,  vapor  shear  is  dominant 
and  the  trends  described   in  this  section  determine   the  heat  transfer.    Figure  IV-6  shows 
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Figure    IV-6.     Effects    of    Vapor    Velocity    on    Mean    Nu    Using    the    Asymptotic 
Interfacial    Shear  Stress  Approximation    Analysis   Method. 
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the  effect  of  F  on  the  mean  heat-transfer    coefficient    for  varying   eccentricity.     As  F  - 
oo,  corresponding    to  free  convection    condensation,    solutions    for  elliptical    and  circular 
tubes   match    those    reported    in  Section    IV. A.     For  the   case  of  a  circular    tube   for  all 
values  of  F,  the  present  numerical   solutions  are  within  0.4%  of  an  equation   by  Rose  [28]: 

&*-*=         M+&728F*  .  (IV-2) 

(1+3.44  Fm  +  F  )1/4 

2.        Two-Phase  Boundary -Layer   Shear  Stress  Approximation 

Vapor  boundary-layer  separation  is  not  predicted  by  the  asymptotic  shear 
stress  approximation  since  the  interfacial  shear  is  based  on  potential  flow  outside  the 
vapor  boundary- layer  which  is  always  positive.  Solution  of  the  two- phase  boundary- 
layer  equations  for  the  condensate  and  vapor,  with  matched  shear  stress  at  the  interface, 
allows  determination  of  the  vapor  boundary-layer  separation  point  and  its  effect  on  the 
condensation  heat  transfer.  The  analysis  also  introduces  a  dimensionless  parameter  G 
which  is  proportional  to  the  condensation  rate.  Figure  IV-5  compares  the  asymptotic 
shear  stress  approximation  with  the  two- phase  boundary  layer  shear  stress  approxima- 
tion for  a  large  condensation  rate  (large  G)  for  a  circular  and  elliptical  tube  (k  =  0.6). 
Over  the  forward  and  middle  sections  of  the  tube,  the  general  trends  are  as  described 
in  Section  IV.B.l  for  the  asymptotic  shear  stress  approximation.  Over  the  rear  of  the 
tube,  however,  as  the  separation  point  is  approached,  the  condensate  film  thickens  more 
rapidly  due  to  the  reduced  vapor  shear  effect  (x  -  0).  It  can  be  seen  that  for  the 
elliptical  tube,  the  separation  point  occurs  further  downstream  than  for  the  circular 
tube  and  thus  this  rapid  thickening  is  delayed.  By  delaying  separation,  more  of  the  tube 
surface  has  a  thinner  film.  This  serves  to  increase  the  heat-transfer  coefficient  for  an 
elliptical    tube  compared   to  a  circular   tube  if  it  were  the  only  effect  considered. 
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In  total,  however,  the  mean  heat-transfer  coefficient  is  affected  by  a 
combination  of  gravity,  vapor  shear  and  boundary- layer  separation.  Figure  IV-7shows 
the  effect  of  F,  G  and  k  on  the  vapor  boundary- layer  separation  point  (xs  ).  As 
previously  explained  (Chapter  III),  solutions  could  not  be  obtained  for  all  possible  F  and 
G  combinations  (particularly  large  F  and  G).  In  general,  as  F  and  G  increase  (less  vapor 
shear  and  higher  condensation  rate),  or  k  decreases  (more  elliptical),  the  vapor 
separation  point  shifts  toward  the  rear  of  the  tube.  For  low  condensation  rates  (small 
G)  corresponding  to  low  vapor  suction,  vapor  separation  occurs  near  the  positions 
obtained  for  single  phase  separation  without  suction.  Little  change  occurs  in  the 
separation  point  location  as  vapor  velocity  is  varied.  For  large  condensation  rates  (large 
G),  the  vapor  boundary-layer  experiences  high  suction  which  shifts  the  separation  point 
downstream  (i.e.,  as  G  increases,  xs  increases).  As  vapor  velocity  is  decreased 
(increasing  F),  the  separation  point  shifts  further  downstream.  As  k  decreases  (more 
elliptical),  the  vapor  boundary- layer  separation  point  moves  further  downstream  when 
compared  to  a  circular  tube  under  all  conditions  of  F  and  G.  At  high  F  and  G,  the  vapor 
boundary- layer  separation  points  c  me  together  at  the  rear  of  the  tube  for  both  the 
elliptical    and  circular    tube. 

Figure  IV-8 shows  the  influence  that  the  above  effects  (vapor  shear,  gravity, 
condensation  rate  and  eccentricity)  have  on  the  mean  heat- transfer  coefficient.  The 
solutions  for  k  -  1  agree  with  the  results  of  Fujii  et  al.  [22]  as  corrected  by  Lee  and  Rose 
[25].  Similar  trends  were  obtained  for  the  elliptical  tubes.  For  large  values  of 
condensation  rate  (large  G),  the  two-phase  boundary-layer  shear  stress  analysis  agrees 
closely  with  the  asymptotic  shear  stress  analysis  method  as  expected  since  the  latter 
assumes  an  infinite  condensation  rate.  For  decreasing  k,  the  results  for  an  elliptical 
tube  tend  to  show  a  reduction    in  heat  transfer    at  low  F  (compared   to  a  circular    tube) 
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VAPOR  SEPARATION  POINT 
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Figure  IV- 7.   Vapor  Separation    Point  Location    with  Varying   F,  G  and  k. 
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the  Two- Phase  Boundary- Layer    Approximation    Analysis    Method. 


and  an  increase  in  heat  transfer  at  high  F,  as  previously  mentioned.  The  cross-over 
point  from  reduction  to  improvement  occurs  at  increasing  values  of  F  as  G  decreases. 
For  large  values  of  F,  the  results  converge  to  the  Nusselt  analysis.  This  shift  in  the 
cross-over  point  is  attributed  to  the  change  in  the  vapor  boundary- layer  separation  point 
with  varying  condensation  and  vapor  velocity  (varying  Gand  F).  As  seen  in  Figure  IV- 
7,  for  large  condensation  rate  (G  =  5.5)  and  low  vapor  velocity  (large  F),  the  large  vapor 
suction  causes  a  delay  in  vapor  boundary- layer  separation  such  that  it  occurs  at 
essentially  the  same  streamwise  location  for  all  eccentricities.  As  velocity  increases  (F 
becomes  smaller)  the  separation  point  shifts  forward  at  different  rates  dependent  on 
eccentricity  which  is  evidenced  by  the  diverging  curves  in  Figure  IV-7c.  Since  early 
vapor  boundary- layer  separation  results  in  a  reduction  in  heat  transfer  due  to  rapid 
thickening    of  the   condensate     film,   the   circular    tube   (k  =  1)  experiences    more   of  a 
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reduction  in  heat  transfer  compared  to  an  elliptical  tube.  The  effect  of  delaying 
separation  is  to  cause  a  thinning  in  the  mean  film  thickness  for  the  tube.  In  the 
balance,  where  the  streamwise  component  of  gravity  and  the  shift  in  vapor  boundary- 
layer  separation  are  positive  effects  (i.e.,  both  contribute  to  a  relative  thinning  of  the 
elliptical  tube  condensate  film  compared  to  a  circular  tube)  and  vapor  shear  is  a 
negative  effect  (i.e.,  lower  shear  contributes  to  a  relative  thickening  of  the  elliptical 
tube  condensate  film  compared  to  a  circular  tube),  delays  in  vapor  separation  move  the 
cross-over  point  from  gravity  dominant  to  shear  dominant  heat  transfer.  For  the  case 
of  low  G  in  Figure  IV-7a,  there  is  no  noticeable  shifting  of  the  vapor  boundary-layer 
separation  point  as  F  is  varied.  Thus  the  vapor  boundary- layer  separation  point 
contribution  does  not  change  with  vapor  velocity  (F)  to  the  balance  and  the  cross-over 
point  occurs  at  lower  velocities  (larger  F).  The  overall  effect  of  vapor  shear  using  this 
analysis  is  typically  a  slight  reduction  (<  2%)  in  the  heat- transfer  coefficient  for  an 
elliptical  tube  compared  to  an  equivalent  surface  area  circular  tube,  dependent  on  the 
magnitude    of  the  condensation    rate  parameter,    G. 

Figures  IV-9and  IV- 10  show  the  effects  of  vapor  boundary-layer  separation 
on  the  local  film  thickness  and  heat- transfer  coefficient  for  a  given  dimensionless 
streamwise  distance.  For  low  F  and  G,  the  separation  point  occurs  relatively  early,  but 
is  significantly  delayed  when  using  an  elliptical  tube.  As  F  increases  (vapor  velocity 
decreases),  the  separation  points  remain  at  their  same  respective  streamwise  locations. 
The  areas  under  the  curves  in  Figure  IV- 10  represent  the  total  heat  transfer.  Bearing 
this  in  mind,  the  comparison  between  elliptical  and  circular  tubes  is  made  a  little  easier. 
At  high  F  (for  given  G),  the  area  under  the  curve  for  an  elliptical  tube  exceeds  that  for 
a  circular  tube  while  at  low  F,  the  opposite  is  true.  However,  for  a  given  eccentricity, 
the  values  of  xs    are  the  same  for  varying  F  (as  noted  earlier). 
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As  G  increases  (increasing  condensation  rate  and  suction),  vapor  boundary-layer 
separation  is  further  delayed  (compared  to  low  G)  and  is  virtually  eliminated  at  high 
F.   Separation    is  always   delayed    for  an  elliptical    tube  compared   to  a  circular    tube. 
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Figure  IV-9  a  and  b.  Dimensionless    Local  Film  Thickness    for  Varying  G,  F  and  k 
from  the  Two- Phase  Boundary    Layer  Shear  Stress  Method  of  Analysis. 
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Figure  IV-9  c  and  d.   Dimensionless    Local  Film  Thickness   for  Varying  G,  F  and  k 
from  the  Two-Phase  Boundary    Layer  Shear  Stress  Method  of  Analysis. 
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Figure  10  a  and  b.  Local  Nu  for  Varying  G,  F  and  k  from  the  Two-Phase  Boundary 
Layer  Shear  Stress  Analysis   Method. 
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Figure   10  c  and  d.   Local  Nu  tor  Varying  G,  F  and  k  from  the  Two- Phase  Boundary 
Layer  Shear  Stress  Analysis    Method. 
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C.        EFFECT  OF  SURFACE    TENSION 

As  noted  in  Section   II. E,  a  condition   exists  for  which  d6*/d<j>  becomes  infinite   at 
some  critical  angle  <{>c.  What  is  actually   happening   can  be  discerned   from  Figure  IV- 11 
which    shows    the   surface    tension    function     f4(<j>)   versus   x     for    varying    eccentricity. 
Surface   tension  causes  a  favorable    pressure  gradient  over  the  front  half  of  the  elliptical 
tube  and  an  adverse   pressure  gradient   over  the  back  half.    The  severity   of  the  pressure 
gradient   is  localized    to  small  regions  at  the  top  and  bottom  of  the  elliptical    tube  where 
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Figure  IV- 11.  Surface   Tension   Function,    f4(4>)  for  Varying   Eccentricity. 

change  in  surface  curvature  is  most  severe.  Therefore,  it  is  most  significant  for  small 
values  of  eccentricity.  The  above  analysis  explains  the  reason  for  condensate  flow 
instability    discussed   in  Section  II. E  and  predicted    by  Equation    (11-90). 


75 


Further  discussion  of  surface  tension  is  restricted  to  those  cases  in  which  a 
solution  could  be  obtained  over  at  least  90%  of  the  tube  surface.  This  restriction  was 
arbitrarily  chosen  based  on  the  small  amount  of  heat  transfer  which  occurs  on  the  lower 
part  of  the  tube  and  would  therefore  have  minimal  impact  on  the  accuracy  of  the 
results.  Figure  IV- 12  shows  the  effects  of  surface  tension  on  the  mean  heat-transfer 
coefficient  for  k  =  0.6  and  0.4  .  The  Bond  number,  Bo,  gives  the  relative  effect  of 
inertia  to  surface  tension.  Values  of  1/Bo  of  0.01  (typical  for  steam)  and  0.001  (typical 
for  refrigerants)  as  well  as  1/Bo  =  0  have  been  shown  on  the  figure.  A  number  of 
interesting  points  can  be  highlighted.  Firstly,  the  surface  tension  effect  is  much  smaller 
for  highly  wetting  refrigerants,  as  expected.  Secondly,  inclusion  of  surface  tension 
leads  to  a  small  decrease  in  the  mean  heat- transfer  coefficient  (<  2%)  over  the  whole 
range  of  F  for  the  practical  range  of  eccentricities,  suggesting  that  any  thinning  of  the 
condensate  film  over  the  top  half  of  the  tube  is  more  than  offset  by  a  thickening  over 
the  lower  half.  As  eccentricity  decreases,  this  discrepancy  is  accentuated.  Finally, 
surface  tension  effects  are  felt  more  in  the  free  convection  region  (high  F)  than  in  the 
forced  convection    region  (low  F)  as  film  thickness   becomes  dominated    by  vapor  shear. 

Krupiczka  [29]  used  the  curvature  of  the  condensate  film  surface  to  analyze 
surface  tension  on  a  horizontal  circular  cylinder.  He  determined  that  the  effect  of 
surface  tension  was  only  significant  over  the  bottom  portion  of  the  tube.  This  fact 
implies  that  since  the  film  is  very  thin  over  the  top  of  the  tube,  its  curvature  is  very 
close  to  that  of  the  tube  surface  itself,  which,  for  a  circular  tube  results  in  no  surface 
tension  effect  (constant  radius  of  curvature).  As  the  condensate  flows  around  the  lower 
half  of  the  circular  tube,  the  film  curvature  no  longer  follows  the  tube  surface.  Rather, 
the  radius  of  curvature  is  increasing  resulting  in  a  small  improvement  in  the  heat- 
transfer   coefficient.     A  similar  analysis    can  be  applied   to  the  elliptical    tube.    Over  the 
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Figure   IV- 12.   Effect  of  Surface    Tension   on  the  Mean  Nu  for  k  =  0.4  and  0.6 
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top  half  of  the  tube  where  the  condensate  film  is  thinnest,  modeling  the  surface  tension 
effect  using  the  curvature  of  the  tube  surface  is  valid.  However,  the  adverse  pressure 
gradient  over  the  lower  half  of  the  tube  is  probably  not  as  severe  as  seen  in  Figure  IV- 
11.  Over  this  portion  of  the  tube,  the  film  is  thicker  and  its  decreasing  surface 
curvature  is  not  as  severe  as  the  tube  surface  curvature.  The  conclusion  drawn  from 
this  analysis  is  that  surface  tension  does  not  have  a  significant  effect  on  the  mean  heat- 
transfer    coefficient    for  an  elliptical    tube. 

D.        EFFECT  OF  PRESSURE    GRADIENT 

The  effect  of  pressure  gradient  was  analyzed  using  the  asymptotic  shear  stress 
approximation  (similar  to  Rose  [28]  for  a  circular  tube).  In  this  analysis,  the  pressure 
gradient  in  the  condensate  film  was  assumed  to  be  due  to  the  pressure  gradient 
impressed  on  the  condensate  by  vapor  potential  flow.  As  in  the  case  of  surface  tension, 
there  is  a  relationship  between  parameters  P  and  F  (Equation  11-89)  for  which  a  solution 
could  not  be  obtained  for  the  entire  tube.  Figure  IV- 13  shows  the  pressure  gradient 
function  versus  x  for  varying  k  which  sets  up  a  favorable  pressure  gradient  over  the 
top  of  the  tube  and  an  adverse  pressure  gradient  over  the  lower  half  of  the  tube.  The 
effect  of  eccentricity  is  to  shift  the  point  of  maximum  pressure  gradient  to  the  front 
and  rear  of  the  tube.  Thus,  <j>c  shifts  downstream  for  an  elliptical  tube  compared  to  a 
circular  tube  for  similar  /'and  F.  This  is  a  result  of  the  more  streamlined  shape  and  the 
smaller   resultant    pressure   drop  over  the  streamwise    length  of  an  elliptical    tube. 

As  in  the  analysis  of  surface  tension,  discussion   is  limited  to  those  cases  for  which 
solutions   could  be  obtained   over  the  entire  tube.   When  pressure  gradient    was  included 
in  the  momentum  equation   of  the  condensate    film,  solutions  agreed  with  the  results  of 
Rose  [28]  for  a  circular    tube,  yielding    a  maximum    5%  decrease    in  the  heat-transfer 
coefficient    over  the  entire  range  of  F  compared   to  the  case  with  pressure  gradient 
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Figure   IV-13.   Pressure    Gradient    Function,    f3(4>),  for  Varying    Eccentricity. 

neglected.  Figure  IV-14  shows  the  mean  heat-transfer  coefficient  for  a  circular  tube 
(k  =  1)  for  varying  F  and  P.  The  favorable  pressure  gradient  increases  the  mean 
condensate  velocity  resulting  in  a  thinner  film  and  an  increased  mean  heat-transfer 
coefficient  over  the  top  of  the  tube.  This  increase  is  offset  by  the  adverse  pressure 
gradient  which  slows  the  condensate  film  velocity  over  the  lower  half  of  the  tube, 
thickens  the  film  and  decreases  the  mean  heat-transfer  coefficient.  For  an  eccentricity 
of  0.6,  Figure  IV- 15  compares  the  effect  of  pressure  gradient  for  an  elliptical  tube  with 
those  for  a  circular  tube  for  varying  F  and  P.  For  both  the  circular  and  elliptical  tubes, 
the  effect  of  pressure  gradient  provides  a  slight  reduction  in  condensation  heat  transfer. 
Figure  IV- 15  indicates  that  the  reduction  in  heat  transfer  was  less  for  the  elliptical  tube 
compared  to  the  circular  tube  as  a  result  of  the  more  favorable  potential  velocity 
distribution    around   the  ellipse. 
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Figure   IV- 14.    Effect   of  Pressure    Gradient   on  the  Mean  Nu  for  a  Circular    Tube 
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Figure    IV- 15.    Effect    of  Pressure    Gradient    on  Mean   Nu,  k  =  0.6,  Relative    to  a 
Circular    Cylinder. 
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E.    EFFECT  OF  VAPOR  PRESSURE  DROP 

The    previous     results     indicate     that     vapor    shear     decreases     the    heat-transfer 
coefficient     for    an    elliptical     tube    when    compared    to   a   circular     tube.      However,    a 
consideration    not  yet  discussed   is  that  the  vapor  pressure   drop  across  an  elliptical    tube 
is    significantly      smaller     than     a    circular     tube    for    the    same    free-stream     velocity. 
Alternatively,     for   a  given    pressure    drop,   the   free-stream    velocity    is  greater    for   the 
elliptical   tube.  Since  very  little  data  is  available    for  two- phase  drag  coefficients   around 
tubes,  single-phase  data  was  used  where  the  drag  coefficient    for  the  tube  is  due  entirely 
to    form  drag.    Values  of  CD  for  a  circular   (k  -  1)  and  an  elliptical    tube  (k  =  0.5)  were 
estimated  to  be  approximately    1.2  and  0.6  respectively   (White  [44]).  For  steam  with  Tsat= 
60°C  condensing   on  a  horizontal    tube  with  Twa|t=  40°C  and  UM=  25  m/s  for  the  circular 
tube,  the  corresponding    vapor   velocity    for  the  elliptical    tube,  for  the  same  pressure 
drop,  was  calculated    to  be  Um=  35  m/s.   The  resulting    values  of  F  are  therefore    0.0257 
and  0.0132  for  the  circular  and  elliptical   tubes  respectively.     Using  the  asymptotic   shear 
stress  assumption    for  both  values   of  F  gives  an  increase    in  the  heat  transfer    for  an 
elliptical     tube    of    16.3%   when    compared    to   a   circular     tube.      Using    the    two-phase 
boundary-layer     shear    stress    assumption     gives    a    corresponding     increase    of    17.1%. 
Therefore,    though   the  effect   of  vapor  shear  alone   results   in  a  reduction    of  the  heat- 
transfer    coefficient    (as  discussed   in  Section  IV.B),  when  the  streamlined    shape  of  the 
elliptical    tube  and  its  effect  on  the  vapor  pressure  drop  is  taken  into  account  (allowing 
for   higher    U.  for   a  given   pressure    drop),   there   is  an   increase    in  the   heat- transfer 
coefficient    when  compared   to  a  circular   cylinder. 
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F.         INSIDE   HEAT  TRANSFER    COEFFICIENT 

So  far   only    the   outside    heat-transfer    coefficient     has   been   discussed.      In   real 

condensers,     it    is    the    overall    heat- transfer     coefficient     which    is    important     and    so 

consideration    must  also  be  given  to  the  single- phase  convective   heat  transfer   occurring 

inside  the  elliptical    tube.  A  preliminary    survey  of  the  literature    (Incropera  and  Dewitt 

[45])  shows  that  for  turbulent   flow  conditions   where  the  Prandtl   number  is  greater  than 

0.5  (which   is  true   for  most  practical    condensers),    correlations    developed    for  circular 

tubes    may    be    used    with    good    accuracy     to    approximate     the    inside    heat- transfer 

coefficient    for   an  elliptical    tube   if  the   hydraulic    diameter    is  used   in  place    of  the 

circular    diameter.     The  hydraulic    diameter    is  given  by 

D    =  t*£  (IV-3) 

P 

where  Ac  is  the  cross-sectional  area  of  the  tube  and  p  is  the  inside  perimeter  of  the  tube. 

The  perimeter   for  an  elliptical    tube  and  an  equivalent   surface  area  circular   tube  is  the 

same.    The  cross-sectional   area  of  an  elliptical    tube  is  given  by 

Ac  =  izab  (IV -4) 

and  is  smaller  than  an  equivalent  surface  area  circular  tube.  Inside  heat  transfer 
correlations  for  turbulent  flow  heat  transfer  inside  a  circular  tube  are  typically  of  the 
form  [45] 
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Thus,  the  inside   heat-transfer    coefficient    is  related    to  the  hydraulic    diameter   by 
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As  an  example,   consider   an  elliptical    tube  with  k  =  0.5  and  a  major  axis  with  a  -  0.01 
m.    The  equivalent    circular    tube  would  have  a  diameter   of  0.0154  m.   The  ratio  of  the 


heat- transfer    coefficient    of  the  elliptical    tube  to  the  circular    tube  is  given  by 
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For  this  simple  analysis,  it  can  be  seen  that  the  inside  heat- transfer  coefficient  of  the 
elliptical  tube  is  3.5%  greater  than  the  equivalent  surface  area  circular  tube.  This 
result,  when  added  to  the  increase  in  the  outside  heat-transfer  coefficient  (for  a  given 
pressure   drop),  would   increase    the  overall   heat- transfer    coefficient    of  the  condenser. 
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V.  CONCLUSIONS 

Analyses  of  laminar  film  condensation  on  a  horizontal  elliptical  tube  were 
conducted  under  conditions  of  free  and  forced  convection,  the  latter  using  both  an 
asymptotic  interfacial  shear  stress  approximation  as  well  as  solution  of  the  two- phase 
boundary- layer  equations.  For  the  asymptotic  shear  stress  approximation,  the  effects 
of  surface  tension  and  pressure  gradient  were  also  included.  Where  possible,  these 
analyses  have  been  validated  against  existing  solutions  for  laminar  film  condensation 
on  a  horizontal    circular    tube  and  a  vertical    flat  plate. 

Whether  or  not  condensation  heat  transfer  is  improved  for  an  elliptical  tube 
compared  to  an  equivalent  surface  area  circular  tube  is  dependent  on  a  balance  of  the 
effects  of  gravity,  vapor  shear  and  vapor  boundary- layer  separation  and  the  influence 
these  have  on  the  thickness  of  the  condensate  film.  Under  quiescent  vapor  conditions 
(no  vapor  shear  or  boundary- layer  separation),  gravity  causes  an  increase  (-10%)  in  the 
heat  transfer  on  an  elliptical  tube  compared  to  a  circular  tube  of  the  same  surface  area 
due  to  an  increase  in  "effective"  gravity.  Under  conditions  of  forced  convection,  both 
shear  stress  approximations  indicate  a  small  decrease  (<  2%)  in  performance  due  to  the 
reduction  in  interfacial  shear  as  a  result  of  the  better  streamlined  shape  of  the  elliptical 
tube  when  compared  to  the  circular  tube  in  the  same  free  stream  vapor  flow.  However, 
when  pressure  drop  effects  are  also  considered,  the  higher  allowable  vapor  velocity  over 
an  elliptical  tube  (for  the  same  pressure  drop  as  a  circular  tube)  results  in  an  increase 
in  the  mean  heat-transfer  coefficient  of  15-20%.  For  conditions  of  mixed  convection, 
the  condensate  film  thickness  is  controlled  by  both  gravity  and  vapor  shear  effects.  The 
cross-over    between    heat    transfer     improvement     (gravity     dominant     flow)    and    heat 


84 


transfer  reduction  (vapor  shear  dominant  flow)  for  the  elliptical  tube  is  dependent  on 
the  vapor  boundary- layer  separation  point.  For  large  condensation  rates  (high  vapor 
suction),  vapor  boundary- layer  separation  is  delayed  and  the  cross-over  point  shifts  to 
higher  velocities  (smaller  F ).  For  low  condensation  rates  (low  vapor  suction),  movement 
of  the  vapor  boundary- layer  separation  point  is  minimal,  resulting  in  a  minor  effect  on 
heat  transfer  and  a  cross-over  point  that  shifts  toward  lower  velocities  (large  F). 
Pressure  gradient  and  surface  tension  each  lead  to  a  small  decrease  in  the  mean  heat- 
transfer    coefficient    (<  2%)  for  an  elliptical    tube. 

In  general,  therefore,  the  outside  heat-transfer  coefficient  of  a  condenser  tube  is 
improved  by  using  an  elliptical  tube  geometry.  Approximation  of  the  single- phase 
inside  heat-transfer  coefficient  using  a  hydraulic  diameter  and  a  Seider-Tate  type 
correlation  indicates  that  the  elliptical  tube  also  has  better  inside  heat  transfer 
performance  compared  to  an  equivalent  surface  area  circular  tube.  Thus,  the  overall 
heat-transfer    coefficient    should  be  improved. 
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VI.    RECOMMENDATIONS 

This  study  represents  an  initial  step  in  determining  if  a  shell-and-tube  condenser 
using  horizontal  elliptical  tubes  has  better  heat  transfer  than  a  comparable  circular  tube 
bundle.  The  effects  of  multi-tube  arrangements  (i.e.,  the  influence  of  adjacent  tubes 
and  the  spacing  between  tubes)  on  vapor  flow  characteristics  needs  to  be  evaluated  to 
determine  if  further  improvements  in  the  overall  heat- transfer  coefficient  can  be 
achieved.  It  is  suspected  that  the  interfacial  shear  of  the  condensate  film  on  elliptical 
tubes  will  be  further  reduced  but  the  corresponding  reduction  in  vapor  pressure  drop 
across  the  condenser  will  result  in  a  net  gain  in  heat  transfer  compared  to  a  circular 
tube  bundle.  An  analysis  similar  to  that  conducted  by  Aoune  and  Burnside  [46,  47]  can 
be  used  to  evaluate    multi-tube  arrangements. 

Further  studies  should  be  conducted  on  a  single  elliptical  tube  to  better  determine 
the  effect  of  surface  tension  on  heat  transfer.  The  surface  tension  model  used  in  this 
study  was  based  only  on  the  curvature  of  the  tube  surface  and  does  not  take  into 
account  the  relatively  thicker  film  over  the  lower  half  of  the  tube.  The  curvature  of 
the  actual  condensate  film  surface  should  be  analyzed  in  a  manner  similar  to  Krupiczka 
[29].  This  analysis,  together  with  the  asymptotic  shear  stress  assumption,  would  involve 
the  solution   of  a  second  order  ODE. 

Analysis  of  the  effect  of  offsetting  the  elliptical  tube  at  an  angle,  a,  with  respect 
to  vapor  flow  and  gravity  must  be  conducted  to  determine  its  effect  on  the  heat- 
transfer  performance  of  the  tube.  This  is  important  since  in  real  condensers,  vapor  will 
approach  the  tube  at  various  angles.  Additionally,  since  the  elliptical  tube  is  a 
streamlined    body,  vapor  flow  in  a  direction   other  than  the  direction   of  the  major  axis 
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may  result  in  a  lifting  force  on  the  tube  which  may  result  in  tube  vibration  and  wear. 
These  effects  would  have  to  be  studied  and  methods  to  minimize  any  resultant  vibration 
must  found.  The  problems  associated  with  joining  elliptical  tubes  to  the  condenser  tube 
sheet  and  baffle    plates  also  needs  to  be  addressed. 

Finally,     experiments     must    be    conducted     on    elliptical     tubes    to    validate     the 
predictions    of  these  theoretical    models. 
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APPENDIX    A  -  DERIVATION    OF  TANGENTIAL    VAPOR   VELOCITY    OVER  AN 
ELLIPTICAL    TUBE  FROM  POTENTIAL    FLOW  THEORY 


The  basic  technique  used  to  determine  fluid  flow  velocity  over  a  body  is  to 
employ  a  conformal  transformation  to  translate  the  body  from  one  complex  plane  into 
another  complex  plane.  From  this,  the  velocity  of  the  transformed  body  can  be  readily 
computed  using  potential  theory.  For  this  particular  application,  the  fluid  velocity 
along  the  surface  of  a  circular  cylinder  is  determined  and  then  transformed  into  the 
complex  plane  in  which  the  cylinder   has  an  elliptical    shape. 

Consider  a  circular  cylinder  in  the  complex  plane,  zl5as  shown  in  Figure  A-l.  The 
surface   of  the  circle  is  defined   by 


z.  =  c  e' 


x?  +  y\  -  c> 


(a) 


(b) 


(A-l) 


The  conformal  transform   to  translate   the  circular  body  from  the  zx  plane  to  an  elliptical 
body  in  the  z2  plane  is  given  by 


h.  =  Zi 


1 


(A -2) 


thus 


h.  =  (*i  +  fy) 


x,  + 


Vi 


1  -  — 
c2 


+  Vi 


1  + 


(A -3) 


The  real  and  imaginary    coordinates    in  the  z2  plane  are  given  by 


X.H--L 


and     y2  =  y, 


l+i- 
c1) 


(A-4) 
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Figure  A-l.  Circular    Cylinder   in  zl  Plane,  Elliptical    Cylinder   in  z2  Plane. 


Substituting    Equation    (A-4)into  Equation    (A-lb)  results   in 


yi 


c 


c   *  — 

c 


=  1 


(A-5) 


which  is  the  expression  for  an  ellipse  in  cartesian  coordinates  where  the  semi- major 
axis,  a,  is  given  by  c  +  (1/c)  and  the  semi-minor  axis,  b,  by  c  -  (1/c).  The  complex 
potential    function    in  the  Zj  plane  for  the  circle  is  given  by 


h<Zj)  =  Um  zxe    l  * 


-<f       Um  c2e'2 


(A-6) 


where  $  is  the  velocity  field  potential   function  and  t|t  is  the  stream  function.    The  fluid 
velocity    in  the  zl  plane  is  determined    from  the  complex  derivative, 


dw 


—  =  Um  e    2 


-«i       U_  cz  e 


i — 

2   „    2 


=   U,    -IV 


1         "1        ' 


(A-7) 
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where  Uj  and  v1  are  the  fluid   velocity   components   in  cartesian    coordinates.     The  fluid 
velocity   in  the  z2  plane  (for  the  ellipse)  is  found  by  the  chain  rule,  thus, 

dw       dwdzx  /a    «\ 

—  = =  u.  -  iv,  (A-o) 

dz,       dzx  dzz 

Combining   Equations    (A-la),(A-2)and  (A-7)  results   in 


dzz  c2e'2e  +  1 


(A-9) 


The  magnitude    of  the  fluid   velocity    is  found  by  multiplying    the  complex  velocity   by 

its  conjugate.    Using  Euler's  formula, 

ei0  =  cos8  +  i  sin6      ,  (A- 10) 

and  with  some  algebraic   manipulation,    the  potential   velocity  of  the  fluid  at  the  surface 
of  the  ellipse  is  given  by 

2  _  dwdw       Ut  c*  (2  +  2cos26)  (A-ll) 


u;  = 


dz^dz^       c*  +  1  +  2c2  cos26 


where   the  bar  denotes   the  complex   conjugate.     Equation    (A-ll)  is  further    simplified 

using  the  trigonometric    identities 

cos20  =  cos26  -  sin26  ,A    ._. 

(A- 12) 

cos26  =  2cos26  -  1      . 

and  the  relationship    between  c,  a  and  b, 

4c2  =  (a+b)2      .  (A"13) 

Thus  the  potential    velocity   is  given  by 

U,  =  U        iUk)  COS^—      .  (A- 14) 

/cos2e  +  fcWe 

The  angle,  8,  is  measured   from  the  xx  (horizontal)    axis.   The  relationship    between  6  and 
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(J>,  the  angle  measured    from  the  v'i  (vertical)    axis,  is  given  by 

6  -■*  -<fr      .  (A- 15) 

2 

Note  that  <J>  is  an  angle  defined  in  the  zx  plane  for  the  circle  and  is  the  parametric   angle 
defined    in  Section   II. A  .   The  potential    velocity    is  thus  given  by 


U    =  U  (!+£)  sin4>         =  v  (1+fc) 


\Jsm2$  +  it2cos24>  <Jl+k*cotfy 


(A- 16) 
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APPENDIX    B  -  COMPUTER    CODE  FOR  ASYMPTOTIC    SHEAR  STRESS   MODE 


*      PROGRAM  MIXED1 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

X 

x  THIS  PROGRAM  NUMERICALLY  SOLVES  THE  PROBLEM  OF  MIXED 

X  CONVECTION  CONDENSATION  WITH  SURFACE  TENSION  ON  AN  ELLIPTICAL 

x  CYLINDER.   THE  PROBLEM  FORMULATION  USES  THE  ASYMPTOTIC  SHEAR 

X  STRESS  APPROXIMATION  ANALYSIS  METHOD. 

x 

x  THE  PROGRAM  USES  AN  AUTHOR  DEVELOPED  FORWARD  STEPPING 

X  ALGORITHM  BASED  ON  A  4TH  ORDER  ADAMS  PREDICTOR-CORRECTOR 

X  METHOD  TO  SOLVE  A  1ST  ORDER  ORDINARY  DIFFERENTIAL  EQUATION. 

x 

x  THIS  PROGRAM  CAN  BE  RUN  WITH  DIMENSIONLESS  PARAMETERS 

X  (F,G,P,  BO)  SPECIFIED  OR  WHERE  THE  PARAMETERS  ARE  CALCULATED 

x  BASED  ON  SPECIFIED  FLUID  PROPERTIES.   MODIFICATIONS  TO  THE 

X  PROGRAM  (REMOVAL  OF  UNNECCESSARY  STEPS  BY  "COMMENTING"  THE 

x  STATEMENT  OUT  OF  THE  PROGRAM)  IS  REQUIRED.    SPECIFIED  FLUID 

x  PROPERTIES  INCLUDE  T  SAT,  T  WALL,  AND  UINF. 

x 

X  NUMERICAL  RESULTS  ARE  DIRECTED  TO  AN  EXTERNAL  DATA  FILE. 

x 

x  MAJOR  VARIABLES  USED  ARE: 

x 

X        THETA      PARAMETRIC  ANGLE 

X        X         STREAMWISE  LENGTH 

X        DEFF       EFFECTIVE  DIAMETER 

x        DEL        CONDENSATE  FILM  THICKNESS 

X        DELND,Z    DIMENSIONLESS  FILM  THICKNESS 

X        F,P,BO,    DIMENSIONLESS  PARAMETERS 

X        JARPRL 

x        UINF      VAPOR  FREE  STREAM  VELOCITY 

X         TS         VAPOR  SATURATION  TEMPERATURE 

X        TW        WALL  SURFACE  TEMPERATURE 

x        DENSL      CONDENSATE  DENSITY 

x        DENSV      VAPOR  DENSITY 

X        CONDL      CONDENSATE  THERMAL  CONDUCTIVITY 

X        VISL       CONDENSATE  VISCOSITY 

X        HFG        LATENT  HEAT  OF  VAPORIZATION 

X        NU        NUSSELT  NUMBER 

X         NUAVG      MEAN  NU  NUMBER 

x        ALPHA      HEAT  TRANSFER  COEFFICIENT 

X         QFLUX      HEAT  FLUX 

x        A         SEMI-MAJOR  AXIS 

X        K         ECCENTRICITY 

x        XFN        PARAMETRIC  RADIUS 

x        G         GRAVITY  FUNCTION 

x        Fl         VELOCITY  FUNCTION 

x        Fll        DERIVATIVE,  VELOCITY  FUNCTION 

x        F2         PRESSURE  GRADIENT  FUNCTION 

x        F3        SURFACE  TENSION  FUNCTION 

x 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

PARAMETERC PI =3. 141592654, GRAV=9. 807) 

REAL  Z( 0 : 180000 ), THETACO: 180000 ),X(0: 180000 ),H, DEFF, XI, X2,X3, 
+NUAVG,NU,NUORE, UINF, TS.TW.TL, DENSL, VISL, TG, CONDL, HFG, B2,DELND, 
+DEL,QFLUX,A,K,RE,SUM,XFN,ALPHA,ALPHAO,DQFLUX,F,PHI, 
+A1.A2,A3,A4,A5,A6,A7,A8.A9,A10,A11,TF,PSAT,TB,TC,TD,R,TE,DENSV, 
+  S I GML  ,  FCN ,  P ,  CI ,  C2 ,  Co  ,  FCO ,  FC1 ,  FC2 ,  FC3 ,  ZP1 ,  ZP2 ,  ZP 3 ,  ZP4  ,  BO , 
+C4,C5,F21,F31,NUCOEFF 

INTEGER  N,J,JJ,IC 

COMMON  /GEOM/ A , K , DEFF/PARAM/F, BO , P 

EXTERNAL  FCN,XFN,G, Gl , Fl , Fll , F2, F3 

0PEN(20,FILE='/MIXED1  OUTPUT  Al • , STATUS= • OLD • ) 
0PEN(30,FILE='/DEL2  DATA  Al ' , STATUS= ■ OLD • ) 
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FILE:  MIXED1    FORTRAN   A 

xxx    INPUT  FLUID  PARAMETERS   XXX 

PRINTX, 'ELLIPSE  OF  THE  FORM  (X/B)xx2+(Y/A)x*2=l ' 

PRINT*, 'INPUT  "A"  DIMENSION  (M)?' 

READ*, A 

PRINT*, 'INPUT  T  SAT  AND  T  WALL  (K)?» 

READ*,TS,TW 

PRINT*, 'INPUT  ANGLE  STEP  SIZE  (DEGREES)?' 

READ*,H 

PRINT*, 'STEP  SIZE,  H  =  «,H 
C     PRINT*, 'INPUT  UINF?' 
C     READ*, UINF 

PRINT*, 'INPUT  K?' 

READ*,K 

PRINT*, 'INPUT  F?' 

READ*,F 

N=180/H+l 
H=PI/180.*H 
BO=0. 
P  =  0. 

***   DETERMINE  FLUID  PROPERTIES    xxxxxxxx 

Al=15. 49217901 
A2=-5. 6783717693 
A3=l. 4597584637 
A4=13. 877000608 
A5=-80. 887673591 
A6=123. 56883468 
A7=-188. 321212064 
A8=660.  91763485 
A9=-1382. 4740091 
A10=1300. 1040184 
All=-449. 39571976 
TF=TS/1000. 

PSAT=1.E+6*EXP(A1+A2/TF+A3*L0G(TF)+TF*A4+A5*TF**2+A6*TF**3+A7* 
+TF**4+A8*TF**5+A9*TF**6+A10*TF**7+A11*TF**8) 
TB=1500./TS 

TC=2.5*L0G(1.-EXP(-TB)) 

TD=.0  015/(1.+.0  001*TS)-.0  0  0942*EXP(TB+TC)*SQRT(1./TB)-.00  04882*TB 
R=461.51 

TE=2.*PSAT/(R*TS) 
DENSV=TE/(1.+SQRT(1.+2.*TD*TE)) 
TL=(2.*TW+TS)/3. 

DENSL=l./(.0012674-TL*(2.02915E-6-TL*3.8333E-9)) 
VISL=2.414E-5*10**(247.8/(TL-140.)) 
TG=TL/273.15 

CONDL  =  -.92407+TG*(2.8395-TG*Q.8007-TG*(.52577-TG*.07344))) 
HFG=3468920.-TS*(57  07.4-TS*(11.5562-TS*.0133103)) 
SIGML=(-.0  003*(TL-273.15)**2-.138*(TL-273.15)+75.6)/1000. 

***   DETERMINE  THETA(I) ,X(I) , DEFF   xxxxxxx* 

XX*   USE  SIMPSON'S  RULE  TO  NUMERICALLY  INTEGRATE  ALONG  SURFACE  xxx 

THETA(0)=0. 
XC0)=0. 

DO  10,  1=1, N 

THETA(I)=THETA(I-1)+H 
X1=A*XFN(THETA(I-D) 
X2=A*XFN((THETA(I-1)+THETA(I))*.5) 
X3  =  A*XFN(THETA(D) 
XCI)=X(I-1)+H/6.*(X1+4.*X2+X3) 
10    CONTINUE 

DEFF=2.*X(N)/PI 

***   DETERMINE  DIMENSIONLESS  PARAMETERS  AS  REQUIRED   xxxx 

C     F=VISLXDEFFXHFG*GRAV/(UINF**2*C0NDL*(TS-TW) ) 
UINF=SQRT(VISL*HFG*GRAV/( F*CONDL*(TS-TW) ) ) 
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FILE:  MIXED1    FORTRAN   A 

C      P=DENSV*HFGXVISL/(DENSl*CONDLX(TS-TW)) 
C      B0=SIGML/(DENSLXGRAVXDEFFXX2) 
RE=DENSLXDEFFXUINF/VISL 

PRINTX, -COMPLETED  THETA, X, DEFF' 

WRITE(20,5000)  »SAT.  TEMP.  (K)  =  ',TS,'WALL  TEMP.  =  »,TW, 
+  »l/BO  =  ',BO,'P  =  »,P,'F  =  »,F,'UINF  (M/S)  =  ',UINF,'REL  =  ',RE 

WRITE(20,5010)  '"A"  LENGTH  CM)  =  ',A,'"K"  ECCENTRICITY  =  • ,K, 
+'DEFF  (M)  =  »,DEFF 

WRITEC20,5020) 

xxx   DETERMINE  INITIAL  CONDITION,  ZO    xxxxxxxx 

B2=SQRT(UINFXDENSL/(DEFFXVISD) 

THETA(0)=l.E-9 

X(0)=l.E-9 

CALL  DRV(F2,0.,F21) 

CALL  DRV(F3,0.,F31) 

Cl  =  (2XFXGK0.)  +  DEFF/AXPXF21+3x(DEFF/A)XX2XFXBOXF31)/3. 

C2=F11(0.) 

C3=-2.XAXXFN(0.)/DEFF 

Z(0)=SQRT((SQRT(C2XX2-4XC1XC3)-C2)/(2.XCD) 

ALPHA0=CONDLXB2/Z(0) 

SUM=0. 

XXX    DETERMINE  DIMENSIONLESS  FILM  THICKNESS,  Z     XXXX 
XXX    ADAMS-MOULTON  PREDICTOR/CORRECTOR  ALGORITHM     XXXX 

FCO=FCN(THETA(0),Z(0)) 

ZP1=Z(0)+HXFC0 

Z(1)=Z(0)+H/2.X(FC0  +  FCN(THETA(1),ZPD) 

FC1  =  FCN(THETA(1),Z(D) 

ZP2=Z(1)+H/2.X(3.XFC1-FC0) 

Z(2)=ZC1)+H/12.X(5XFCN(THETA(2),ZP2)+8.XFC1-FC0) 

FC2=FCN(THETA(2),ZC2)) 

ZP3=ZC2)+H/12.X(23.XFC2-16.XFC1+5.XFC0) 

ZC3)=ZC2)+H/24.x(9.xFCN(THETAC3),ZP3)+19.xFC2-5.XFCl+FC0) 

FC3=FCN(THETA(3),Z(3)) 

DO   50, 1=4, N 

ZP4=Z(  I-D+H/24  .  X(55  .  XFC3-59  .  XFC2+37  .  XFC1-9  .  XFCO ) 
Z(I)=Z(I-l)+H/720.x(251.XFCN(THETA(I),ZP4)+646.xFC3-264.x 
+        FC2+106.XFC1-19.XFC0) 
FC0=FC1 
FC1=FC2 
FC2=FC3 

FC3  =  FCN(THETA(I),Z(D) 

C4=(2.XFXG(THETA(I))+DEFF/AXPXF2(THETA(I))+3.X(DEFF/A)XX2XF 
+        XB0XF3(THETA(I)))XZ(I)XX3+F1(THETA(I))XZ(I) 
C5=C4+F1(THETA(I))XZ(I) 
IFCC4.LE.0.)    THEN 

PRINTX,  'C4  LIMIT' 
IC  =  I 

GO  TO  60 
ENDIF 
IFCC5.LE.0.)  THEN 

PRINTX, 'C5  LIMIT' 
IC=I 

GO  TO  60 
ENDIF 
50     CONTINUE 
IC=N-1 

60     PRINTX, 'COMPLETED  Z  DETERMINATION' 

xxx    DETERMINE  HEAT  TRANSFER  PROPERTIES    xxxxxxxxx 

J  =  0 

JJ  =  0 

DO  100, 1=1, IC 
DELND=Z(I) 
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FILE:  MIXED1    FORTRAN   A 

DEL=DELND/B2 

ALPHA=CONDL/DEL 

SUM=SUM+CALPHA0XAXXFNCTHETACI-1))+ALPHAXAXXFNCTHETACI)))XHX.5 

ALPHAO=ALPHA 

IF(JJ.EQ.J)  THEN 

QFLUX=ALPHAX(TS-TW) 

NU=DEFF/DEL 

DQFLUX=1./DELND 

HRITEC20,5030)  THETA C  I)*180  ./PI,  X(  I)X2  ./(  DEFFXPI ) , 
+         DELX10  0  0.,QFLUX/10  00. , DELND, DQFLUX, NU 

WRITEC30,x)  X(I)X2./(DEFF*PI),DELND 

J=J+1./180.XN+1 
ENDIF 
JJ=I+1 
100    CONTINUE 

WRITEC20,5030)  THETAC IOX180  ./PI  ,XC  IC)X2  ./C  DEFFXPI ),  DELxl 000  .  , 
+QFLUX/ 1000., DELND, DQFLUX, DEFF/DEL 
WRITEC30,X)  XCIC)X2./(DEFF*PI), DELND 

NUAVG=SUMX2/(PIXC0NDL) 

WRITEC20,5040)  'NU/SQRTCRE)  =  ' , NUAVG/SQRTCRE) , 'NUAVG  =  ',NUAVG 

5000   F0RMATC1X,T5,70C'X»)  //  IX, T5, 'HEAT  TRANSFER  PROPERTIES'  //  2(1X, 

+T10,A17,F6.2  /  ),4C1X,T10,A17,E9.3  /)) 
5010   FORMATC  /  IX, T5, • ELLIPSE  GEOMETRY'  //  3( 1X,T10 , A19,  F6 . 4  /  )) 
5020   FORMATC  /  IX, T7, 'THETA ' , T15, ' DIM  X' , T25, » DEL ' , T32, 'HT  FLUX', 

+T41,'DIM  DEL',T52, *DIM',T58, 'NUSS  #*  /  T7, ' CDEG) « , T25, • CMM) » , T32, 

+'CKW/M2)',T50, *HT  FLUX'  /  T6 , 7C » = • ) , T15,6 C ' = • ) , T23, 7C '=' ) , T32, 

+7C'=,),T41,7C'='),T50,7C'=I),T58,6C '=')  /) 
5030   FORMATC1X,T6,F7.3,T15,F6.4,T23,F7.4,T32,E8.2,T41,F7.4,T50,F7.4, 

+T58,E9.3  ) 
5040   FORMATC  //  3C IX, T5, A17, F8 . 4  /)) 
6000   FORMATC1X,F4.1,F10.3,F7.3) 

END 
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

FUNCTION  XFNCPHI) 
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

X 

x     THIS  FUNCTION  SUBPROGRAM  IS  USED  TO  DETERMINE  STREAMWISE  LENGTH 

x     AS  A  FUNCTION  OF  THETA.   XFN  REPRESENTS  PARAMETRIC  RADIUS. 

X 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

REAL  XFN,PHI,A,K,DEFF 
COMMON  /GEOM/A,K,DEFF 

XFN=SQRTC CKXX2+. 25XC CKXX2-1 . )XSINC2XPHI) )XX2)/C C COSC PHI) )XX2 
++CKXSINCPHI))XX2)) 
END 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

FUNCTION  GCPHI) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

X 

x     THIS  FUNCTION  SUBPROGRAM  REPRESENTS  THE  GRAVITY  FUNCTION. 

X 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

REAL  G,PHI,A,K,DEFF 
COMMON  /GEOM/A,K,DEFF 

G=SINCPHI)/SQRTCCSINCPHI))XX2+CKXC0SCPHI))XX2) 
END 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
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FILE:  MIXED1    FORTRAN   A 

FUNCTION  G1CPHI) 
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

X 

x      THIS  FUNCTION  SUBPROGRAM  REPRESENTS  THE  STREAMWISE  DERIVATIVE 

X      OF  THE  GRAVITY  FUNCTION. 

x 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

REAL  G1,PHI,A,K,DEFF 

COMMON/GEOM/A, K, DEFF 

Gl=Kxx2xC0S(PHI)/SQRT(CCSIN(PHI))xx2+(KXC0S(PHI))xx2)xx3) 

END 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
FUNCTION    FKPHI) 

*XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

X      THIS  FUNCTION  SUBPROGRAM  REPRESENTS  THE  VAPOR  VELOCITY 
X      FUNCTION 

*xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

REAL  F1,PHI,A,K,DEFF 

COMMON  /GEOM/A,K,DEFF 

F1  =  (1.+K)XSIN(PHIVSQRTC(SIN(PHI))XX2+(KXC0SCPHI))XX2) 

END 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

FUNCTION  Fll(PHI) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

X 

X      THIS  FUNCTION  SUBPROGRAM  REPRESENTS  THE  STREAMWISE  DERIVATIVE 
X      OF  THE  VAPOR  VELOCITY  FUNCTION, 
x 
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

REAL  F11,PHI,A,K,DEFF 

COMMON/GEOM/A, K, DEFF 

F11=(1.+K)XKXX2XC0S(PHI)/SQRT(((SIN(PHI))XX2+(KXC0S(PHI))XX2)XX3) 

END 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

FUNCTION  FZ(PHI) 
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

X 

x      THIS  FUNCTION  SUBPROGRAM  REPRESENTS  THE  PRESSURE  GRADIENT 

X      FUNCTION. 

X 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

REAL  F2,PHI,A,K,DEFF 
EXTERNAL  XFN 
COMMON/GEOM/A, K,  DEFF 

F2=((1.+K)XK)XX2XSIN(2XPHI)/((SIN(PHI))XX2+(KXC0S(PHI))XX2) 
+/XFN(PHI) 
END 

*xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

FUNCTION  F3CPHI) 

KXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

X 

X      THIS  FUNCTION  SUBPROGRAM  REPRESENTS  THE  SURFACE  FUNCTION. 

X 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
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FILE:  MIXED1    FORTRAN 


REAL  F3,PHI,A,K,DEFF 
EXTERNAL  XFN 
COMMON/GEOM/A,K,DEFF 

F3=KX(1.-KXX2)XSIN(2XPHIVSQRT(C(SIN(PHI))XX2+(KXC0SCPHI))XX2) 
+XX5)/XFNCPHI) 
END 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
FUNCTION  FCN(PHI,Z) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

X 

X      THIS  FUNCTION  SUBPROGRAM  IS  USED  IN  A  PROBLEM-SOLVER 

X      SUBROUTINE  TO  SOLVE  THE  O.D.E. 

X  DZ/DX  =  F(PHI,Z) 

X 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

REAL  F21,F31,PHI,FCN1,FCN2,FCN,A,K,DEFF,F,B0,P 

COMMON  /GEOM/A, K, DEFF/PARAM/F, BO, P 

EXTERNAL  XFN,G,G1, Fl, Fll, F2, F3 

CALL    DRV(F2,PHI,F21) 
CALL    DRVCF3,PHI,F31) 

FCN1=(2.XFXGCPHI)+DEFF/AXPXF2CPHI)+3.X(DEFF/A)*X2*F*B0XF3CPHI)) 
+XZXX3+FKPHDXZ 

FCN2=2.XAXXFN(PHI)/DEFF-(2.XFXG1(PHI)+DEFF/AXPXF21+3.X(DEFF/A)XX2 

+XFXB0XF31)XZXXV3.-F11(PHI)XZXX2 

FCN=FCN2/FCN1 

END 
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

SUBROUTINE   DRV(FN, PHI, FN1) 
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

X 

X     THIS  SUBROUTINE  NUMERICALLY  DIFFERENTIATES  A  FUNCTION. 

X 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

REAL  FN,PHI,FN1,EPS 

FN1=( FN( PHI+1 . E-6 )-FN(PHI-l . E-6 ) )/2 . E-6 

RETURN 

END 
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APPENDIX    C  -COMPUTER   CODE  FOR  TWO-PHASE  BOUNDARY    LAYER 
MODEL 

x  PROGRAM  VAPBL1 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

X 

X  THIS  PROGRAM  NUMERICALLY  SOLVES  THE  PROBLEM  OF  MIXED 

X  CONVECTION  CONDENSATION  WITH  SURFACE  TENSION  ON  AN  ELLIPTICAL 

X  CYLINDER.   THE  PROBLEM  FORMULATION  USES  THE  SOLUTION  OF  THE 

X  TWO-PHASE  BOUNDARY  LAYER  EQUATIONS  FOR  THE  CONDENSATE  AND  VAPOR. 

X 

x  THE  PROGRAM  USES  AN  IMSL  NUMERICAL  PROBLEM  SOLVER,  DIVPAG, 

x  TO  SOLVE  A  1ST  ORDER  ORDINARY  DIFFERENTIAL  EQUATION.   THE 

x  PROBLEM  SOLVER  IS  WRITTEN  IN  DOUBLE  PRECISION  WHILE  THE  MAIN 

X  PROGRAM  IS  IN  SINGLE  PRECISION.   THE  MAIN  PROGRAM  MUST  BE 

x  COMPILED  IN  DOUBLE  PRECISION  (AUTO  DOUBLE). 

x 

x  THIS  PROGRAM  CAN  BE  RUN  WITH  DIMENSIONLESS  PARAMETERS 

X  (F,G)  SPECIFIED  OR  WHERE  THE  PARAMETERS  ARE  CALCULATED 

x  BASED  ON  SPECIFIED  FLUID  PROPERTIES.   MODIFICATIONS  TO  THE 

X  PROGRAM  (REMOVAL  OF  UNNECCESSARY  STEPS  BY  "COMMENTING"  THE 

x  STATEMENT  OUT  OF  THE  PROGRAM)  IS  REQUIRED.    SPECIFIED  FLUID 

x  PROPERTIES  INCLUDE  T  SAT,  T  WALL,  AND  UINF. 

x 

X  IF  THE  ALGORITHM  HAS  DIFFICULTY  LOCATING  THE  VAPOR  SEPARATION 

X  POINT,  THEN  THE  PROGRAM  NEEDS  TO  BE  MODIFIED  TO  INPUT  THE 

x  PARAMETRIC  AND,  PHI,  FOR  THE  SEPARATION  POINT. 

x 

X  NUMERICAL  RESULTS  ARE  DIRECTED  TO  AN  EXTERNAL  DATA  FILE. 

x 

x  MAJOR  VARIABLES  USED  ARE: 

x 

x        PHI        PARAMETRIC  ANGLE 

x        X         STREAMWISE  LENGTH 

x        DEFF       EFFECTIVE  DIAMETER 

x        DEL        CONDENSATE  FILM  THICKNESS 

X        DELND      DIMENSIONLESS  FILM  THICKNESS 

x        2         DIMENSIONLESS  VAPOR  MOMENTUM  THICKNESS 

x        F.JARPRL   DIMENSIONLESS  PARAMETERS 

x        UINF      VAPOR  FREE  STREAM  VELOCITY 

X         TS         VAPOR  SATURATION  TEMPERATURE 

X        TW        WALL  SURFACE  TEMPERATURE 

X        DENSL      CONDENSATE  DENSITY 

x        DENSV      VAPOR  DENSITY 

X        CONDL      CONDENSATE  THERMAL  CONDUCTIVITY 

X        VISL       CONDENSATE  VISCOSITY 

X        VISV       VISCOSITY  OF  VAPOR 

X        HFG        LATENT  HEAT  OF  VAPORIZATION 

X         NU         NUSSELT  NUMBER 

X        REV        VAPOR  REYNOLDS  NUMBER 

X        RE        TWO-PHASE  REYNOLDS  NUMBER 

X         NUAVG      MEAN  NU  NUMBER 

x        ALPHA      HEAT  TRANSFER  COEFFICIENT 

X        QFLUX      HEAT  FLUX 

X         A  SEMI-MAJOR  AXIS 

x        K  ECCENTRICITY 

X         XFN        PARAMETRIC  RADIUS 

X         G  GRAVITY  FUNCTION 

x        Fl         VELOCITY  FUNCTION 

X         Fll        DERIVATIVE,  VELOCITY  FUNCTION 

X        F2        PRESSURE  GRADIENT  FUNCTION 

X         F3         SURFACE  TENSION  FUNCTION 

x        TAU        INTERFACIAL  SHEAR 

X         KAPPA      PRESSURE  GRADIENT  FORM  PARAMETER 

X         KAPPAA     SUCTION  FORM  PARAMETER 

x 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXKXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

PARAMETER( PI =  3. 141592654, GRAV=9. 807, NEQ=2, NEQ1 =1 , NPARAM=50 ) 

REAL    PHKO:  180000),  DEL  (0:  180000), X(0:  180000), A,  TS.TW, 
+H, A1,A2, A3, A4,A5,A6,A7,A8,A9,A10,A11,TF,PSAT,TB,TC,TD,R,TE, DENSV, 
+TL, DENSL, VI SL, VI SV,TG, CONDL, HFG, SIGML,F,K, XI, X2,X5, DEFF, UINF, 
+REV,RE,JARPRL,DEL0O,DEL0,Z0,TAU,DTAU,FDEL0,DFDELO,DELl,Zl,ALPHAO, 
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FILE:  VAPBL     FORTRAN   A 

+  SUM,AACl,n,FCN,FCNJ,HINIT,PARAM(NPARAM),KAPPA,KAPPAA,DKAPPA, 
+TOL,XEND,X3EG,YC2),DDEL,DZ,Y1(1),NU,NUAVG,NUORE,DKAPPAA,XO,DXO, 
+DELND, DELL, ALPHA, QFLUX, DQFLUX, D2UPHI , DUPHI, UPHI 

INTEGER  N,JJ,II,IDO,IMETH 

COMMON  /GEOM/A,K,DEFF/PARAM/JARPRL,F 

EXTERNAL  ZFN, DELFN, XFN, G,G1, Fl, Fll , NUSSELT, DRV, FCN, FCNJ , IVPAG, 
+SSET 

OPEN(20,FILE='/VAPOR1  OUTPUT  Al ' , STATUS= 'OLD' ) 
OPEN(30,FILE=WAPOR2  OUTPUT  Al ' ,  STATUS= 'OLD' ) 

PRINT*, 'ELLIPSE  OF  THE  FORM  (X/B)SX2+(Y/A)XX2=1 ■ 

PRINTS, 'INPUT  "A"  DIMENSION  (M)?' 

READS, A 

PRINTX, 'INPUT  T  SAT  AND  T  WALL  (K)?» 

READS, TS,TW 

PRINTX, 'INPUT  ANGLE  STEP  SIZE  (DEGREES)?' 

READS, H 

PRINT*, 'STEP  SIZE,  H  =  ',H 
C     PRINTS, 'INPUT  PHI  AT  SEPARATION?' 
C     READ*, PHIS 

PRINTS, 'INPUT  UINF?' 

READS, UINF 

PRINTS, 'INPUT  K?' 

READ*,K 
C      PRINT*, 'INPUT  G?' 
C      READX,JARPRL 

N=180/H+l 
PRINTS, 'N-1=',N-1 
H=PI/180.SH 

XXX   DETERMINE  FLUID  PROPERTIES   XXXXXXXX 

Al  =  15. 49217901 
A2=-5. 6783717693 
A3=l. 4597584637 
A4=13. 877000608 
A5=-80. 887673591 
A6=123. 56883468 
A7=-188. 321212064 
A8=660. 91763485 
A9=-1382. 4740091 
A10=1300. 1040184 
All=-449. 39571976 
TF=TS/1000. 

PSAT=1.E+6XEXP(A1+A2/TF+A3XL0G(TF)+TFSA4+A5STFSS2+A6XTFSS3+A7S 
+TFXX4+A8XTFXX5+A9XTFXX6+A10XTFXX7+A11XTFSS8) 
TB=1500./TS 

TC=2.5SL0G(1.-EXP(-TB)) 

TD=.0  015/(1.+.0  0  01STS)-.000942SEXP(TB+TC)SSQRT(1./TB)-.0  0  04882STB 
R=461.51 

TE=2.*PSAT/(RXTS) 
DENSV=TE/(1.+SQRT(1.+2.STDSTE)) 
TL=(2.*TW+TS)/3. 

DENSL=l./(.0012674-TLX(2.02915E-6-TLX3.8333E-9)) 
VISL=2.414E-5X10XX(247.8/(TL-140.)) 
VISV=-4.478415E-6+TSX(5.0216E-8-1.579E-ll*TS) 
TG=TL/273.15 

CONDL=-.92407+TGS(2.8395-TGS(1.8007-TGS(.52577-TGS.07344))) 
HFG=3468920.-TSX(57  07.4-TSX(11.5562-TSX.0133103)) 
SIGML=(-. 0003S(TL-273.15)SS2-. 138S( TL-273. 153+75. 6 )/l 000. 

xxx   DETERMINE  PHICI) , X(I) , DEFF   xxxxxxxx 

xxx   USES  SIMPSON'S  RULE  TO  NUMERICALLY  INTEGRATE  ALONG  SURFACE   *sx 

PHI(0)=0. 
X(0)=0. 
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FILE:  VAPBL     FORTRAN   A 

DO  10,  1=1, N 

PHI(I)=PHI(I-1)+H 
X1=XFN(PHI(I-1)) 
X2=XFN((PHI(I-1)+PHI(I))X.5) 
X3=XFN(PHI(I)) 

X(I)=X(I-1)+H/6.X(X1+4.XX2+X3) 
10     CONTINUE 

DEFF=2.XXCN)/PI 

PRINT*, 'COMPLETED  DE  DETERMINATION,  DE=',DEFF 

xxxx     DIMENSIONLESS  PARAMETER  DETERMINATION     xxxx 

F=VISLXDEFFXHFG/(C0NDLxUINFXX2X(TS-TW)) 

REV=UINF*DENSV*DEFF/VISV 

RE=UINFXDENSLXDEFF/VISL 

JARPRL=CONDL*CTS-TW)/HFGXSQRTCDENSL/(DENSVXVISLXVISV)) 

PRINT*, 'COMPLETED  PARAMETER  DETERMINATION* 

WRITE(20,5000)  'SAT.  TEMP.  (K)  =  ',TS,*WALL  TEMP.  =  «,TW, 
+»F  =  ',F, 'UINF  (M/S)  =  ',UINF,«G  =  »,JARPRL 

NRITE(20,5010)  '"A"  LENGTH  (M)  =  ',A,'"K"    ECCENTRICITY  =  ',K, 
+'DEFF  CM)  =  ',DEFF 

XXX  DETERMINE  INITIAL  VALUES  OF  Z  AND  DEL  *** 

x*x  USES  INTIAL  VALUE  OF  DEL  FROM  SHEKRILADZE-GOMELAURI  xx* 

xxx  FOR  MIXED  CONVECTION  AS  STARTING  POINT  TO  FIND  Z  THEN  xxx 

xxx  REDETERMINES  DEL  FOR  THE  PREDICTED  VALUE  OF  Z  USING  xxx 

XXX  NEWTON'S  METHOD,  CHECKS  FOR  CONVERGENCE  AND  ITERATES  xxx 

DEL00=SQRT( . 75XK/FX(SQRT( ( ( 1 .+K)/K)XX2+16 . XFXA/(3 . XDEFF) ) 
+-C1.+KVK)) 

DEL0=DEL00 

Z0=(-.07  35XXFN(0.)XJARPRL/DEFF/F11(0.)/DEL0+SQRT((.0735X 
+XFNCO. )XJARPRL/DEFF/F11(0.)/DEL0)XX2+.0735XXFN(0.)/DEFF/ 
+F11(0.)))XX2 

UPHI=F1(0.) 

DUPHI=F11(0.) 

CALL    DRV(F11,0.,D2UPHI) 

XO=XFN(0.) 

CALL    DRV(XFN,0.,DXO) 

20     KAPPA=DEFF/X0XDUPHIXZ0 

30     KAPPAA=.0682+.17  4XJARPRLXSQRTCZ0)/DEL0 

DKAPPA=DEFFXZ0XD2UPHI/X0-DEFFXDUPHIXZ0XDX0/X0XX2 
DTAU=6.44XSORT(KAPPAXKAPPAA+KAPPAAXX2)XDUPHI/SQRT(ZO)+3.22XUPHI 
+XKAPPAAXDKAPPA/SQRT(Z0)/SQRT(KAPPAXKAPPAA+KAPPAAXX2) 

FDEL0  =  X0/DEFF-F/3.XG1(0.  )xDEL0XX4-DTAUXDEL0xx3/(  <+.  XJARPRL  ) 
DFDEL0=-4.XFXG1(0. )xDEL0XX3/3 . -3 . XDTAUXDEL0XX2/C4 . XJARPRL ) 
DEL1=DEL0-FDEL0/DFDEL0 
IF((ABS(DELl-DEL0).GT.(.001XDEL0)).AND.(FDEL0.GT.l.E-7))  THEN 

DEL0=DEL1 

GO  TO  30 
ENDIF 

Z1=(-.0735XXFN(0. ) XJARPRL/ ( DEFFXF11 ( 0 . )XDEL1 )+SQRT( ( .0735* 
+XFNC0. )XJARPRL/(DEFFXF11(0. )XDEL1 ) )XX2+ . 07  35XXFN( 0 . )/(DEFFX 
+F11C0. ))))*X2 

IF  ((ABSCDELl-DEL00).GT.(.001XDEL0O)).OR.(ABS(Zl-Z0).GT.( .OOlx 
+Z0)))  THEN 

Z0  =  Z1 

DEL00=DEL1 

DEL0=DEL1 

GO  TO  20 
ENDIF 

DEL(0)=DEL1 
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FILE:  VAPBL     FORTRAN 


ALPHA0=1./DEL(0) 
SUM=0. 

PRINT*, 'COMPLETED  Z(0),DELC0)  DETERMINATION' , »    Z=»,Z1, 
+'DEL=',DEL(0) 

xx*    EVALUATION  OF  Z(I)  AND  DEL(I)  **** 

XX*    USES  IMSL  ROUTINE  DIVPAG  xxxx 

HINIT=l.E-<+ 
IMETH=2 

CALL  SSET(NPARAM,0.0,PARAM,1) 
PARAM(1)=HINIT 
PARAM(12)=IMETH 
XBEG=PHICO) 
YC1)=Z1 
YC2)=DEL(0) 
T0L=l.E-3 
IDO=l 

DO  100,I=1,N-1 
C        IF(PHKI).GT.PHIS)  GO  TO  190 

XEND=PHI(I) 

CALL  DIVPAG(IDO,NEQ,FCN,FCNJ,AA,XBEG,XEND,TOL,PARAM,Y) 

DEL(I)=Y(2) 

KAPPA  =  DEFF/XFN(PHI(I))*F11(PHKI))*Y(1) 

KAPPAA=.0682+.17  4*JARPRL*SQRT(Y(1))/Y(2) 

TAU  =  6.<+4xFl(PHICI))*SQRT(KAPPAXKAPPAA+KAPPAA**2)/SQRT(Y(l)) 

WRITEC30,*)  'PHI=',PHICI)X180/PI, 'TAU=',TAU 

11  =  1 

X*X    CHECK  FOR  VAPOR  BOUNDARY  LAYER  SEPARATION      xxxx 

IFCTAU.LE.O.)  GO  TO  190 
100    CONTINUE 

190    PRINT*, 'SEPARATION  ANGLE  AT  PHI  =  ' ,  PHK  II)*180  ./PI 
PRINT*, 'TAU=',TAU, ' DEL= ' , DELCII) 

ID0  =  3 

CALL  DIVPAG(IDO,NEQ,FCN,FCNJ,AA,XBEG,XEND,TOL,PARAM,Y) 

***  USES  A  SIMPLE  NUSSELT  ANALYSIS  TO  COMPUTE  FILM  THICKNESS  **** 
xxx  DOWNSTREAM  OF  THE  SEPARATION  POINT.  PROBLEM  SOLVER  IS  xxxx 
X**    IMSL  DIVPAG  **** 

IFCII.NE.(N-l))  THEN 
XBEG=PHI(II) 
ID0  =  1 

Y1(1)=DEL(II) 
DO  200,  1=11+1, N-l 
XEND=PHI(I) 

CALL  DIVPAG(ID0,NEQ1, NUSSELT, FCN J, AA,XBEG, XEND, TOL,PARAM 
+         ,Y1) 

DELCI)=Y1(1) 
200      CONTINUE 
ENDIF 
ID0  =  3 
CALL  DIVPAG(ID0,NEQ1, NUSSELT, FCNJ, AA,XBEG, XEND, TOL,PARAM, YD 

PRINT*, 'COMPLETED  ZCI),  DELCI)  DETERMINATION* 
WRITE(20,5050)  PHI(II)*180 ./PI 
WRITE(20,5020) 

***    DETERMINE  HEAT  TRANSFER  PROPERTIES    ********* 

JJ  =  0 

DO  300, 1=1, N-l 

DELND=DEL(I) 

DELL=DELND*SQRT(VISL*DEFF/(DENSL*UINF)) 

ALPHA=1./DEL(I) 

SUM=SUM+(ALPHA0*XFN(PHI(I-1))+ALPHA*XFN(PHI(I)))*H*.5 

ALPHA0=ALPHA 
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FILE:  VAPBL     FORTRAN   A 

IF(CJJ.EQ.I)  .OR.UJ.EQ.O))  THEN 
QFLUX=AIPHAX(TS-TW) 
NU=DEFF/DELL 
DQFLUX=1./DELND 
JJ=JJ+10./130.*N+1 
WRITE(20,50  30)  PHI ( I JX18  0/PI , X( I )X2/(DEFF*PI) , DELL 
+         *1000,QFLUX/1000,DELND,DQFLUX,NU 
ENDIF 
300   CONTINUE 

WRITEC20,5030)  PHKN-1  )X18  0  ./PI ,  XCN-1  )x2./(  DEFFXPI) ,  DELLX1000  . , 
+QFLUX/1000.,DELND,DQFLUX,DEFF/DELL 
WRITE(30,X)  XCN-1)X2./(DEFFXPI),DELND 

NUAVG=SUMX2/(PI*DEFF) 

WRITEC20,5040)  «NU/SQRT(RE)  =  • , NUAVG/SQRTCRE) 

5000   FORMAT(lX,T5,70('x»)  //  IX, T5, 'HEAT  TRANSFER  PROPERTIES'  //  2C1X, 

+T10,A17,F6.2  /  ),5(1X,T10,A17,E9.3  /)) 
5010   FORMATC  /  IX, T5,  '  ELLIPSE  GEOMETRY'  //  3(  1X,T10,  A19  ,  F6  .  4  /  )) 
5020   FORMATC  /  IX,  T7,  'THETA • ,  T15,  ' DIM  X' , T25, ' DEL ' , T32, 'HT  FLUX', 

+T41,'DIM  DEL',T52, 'DIM',T58, *NUSS  #•  /  T7, ' ( DEG) » , T25, * (MM) • , T32, 

+'(KW/M2)',T50, 'HT  FLUX'  /  T6,7( ' = • ) , T15,6C ' = * ) , T23, 7( • = ' ) , T32, 

+7C'='),T41,7('=*),T50,7C'='),T58,6(I=')  /) 
5030   F0RMAT(1X,T6,F7.3,T15,F6.4,T23,F7.4,T32,E8.2,T41,F7.4,T50,F7.4, 

+T58,F9.1) 
5040   FORMATC  //  3C IX, T5, A17, F8 . 4  /)) 
5050   F0RMATC1X,T6, 'VAPOR  SEPARATION  AT  PHI  =  «,F6.2) 

END 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

FUNCTION  XFNCPHI) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

X 

X     THIS  FUNCTION  SUBPROGRAM  IS  USED  TO  DETERMINE  STREAMWISE  LENGTH 

X     AS  A  FUNCTION  OF  PHI. 

X 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

REAL  XFN,PHI,A,K,DEFF 
COMMON  /GEOM/A,K,DEFF 
XFN= AXSQRT C C KXX2+. 25XCCKXX2-1.)XSINC2XPHI))X*2)/CC COSC PHI ))XX2 

++(K*SIN(PHI))*X2)) 
END 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

FUNCTION  G(PHI) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
X 

X     GRAVITY  FUNCTION 
x 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

REAL  G,PHI,A,K,DEFF 
COMMON  /GEOM/A,K,DEFF 

G=SIN(PHI)/SQRT((SIN(PHI))XX2+(KXC0S(PHI))XX2) 
END 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

FUNCTION   GKPHI) 
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

X 

X  STREAMWISE   DERIVATIVE   OF   GRAVITY    FUNCTION 

X 
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FILE:  VAP3L     FORTRAN   A 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

REAL    G1,PHI,A,K,DEFF 

COMMON/GEOM/A,K,DEFF 

G1=KXX2XC0S(PHI)/SQRT(((SIN(PHI))XX2+(K*C0S(PHI))XX2)XX3) 

END 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

FUNCTION    FKPHI) 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

X      VAPOR  VELOCITY  FUNCTION 

x 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

REAL  F1,PHI,A,K,DEFF 
COMMON  /GEOM/A,K,DEFF 

F1=(1.+!OXSIN(PHI)/SQRTCCSINCPHI))XX2+(KXCOSCPHI))XX2) 
END 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

FUNCTION  Fll(PHI) 
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

X 

X      STREAMWISE  DERIVATIVE  OF  VELOCITY  FUNCTION 

x 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

REAL  F11,PHI,A,K,DEFF 

COMMON/GEOM/A,K,DEFF 

F11=(1.+K)XKXX2XC0S(PHI)/SQRT(((SIN(PHI))XX2+(KXC0S(PHI))XX2)XX3) 

END 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

FUNCTION  DELFN(PHI,Z,DEL,DZ) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

X 

*     STREAMWISE  DERIVATIVE  OF  CONDENSATE  FILM  THICKNESS 

x 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

REAL  DELFN,PHI,Z,DEL,DZ,DELFN1,DELFN2,DELFN3,X,FF1,FF11,FF12,DXFN 
+,A,K,DEFF,JARPRL,F,KAPPA,KAPPAA,XX,FKAPPA 

EXTERNAL  XFN, Fl , F11,G,G1 

COMMON/GEOM/A,K,DEFF/PARAM/JARPRL,F 

IFCPHI.EQ.O.D+0)  THEN 

DELFN=0.D+0 
ELSE 

FF1=F1(PHI) 
FF11=F11(PHI) 
CALL  DRV(F11,PHI,FF12) 
X=XFN(PHI) 

CALL    DRV(XFN,PHI,DXFN) 
KAPPA=DEFF*FF11*Z/X 

KAPPAA=  .  0682+  .  174XJARPRLXSQRTCZVDEL 
XX=KAPPAAXKAPPA+KAPPAAXX2 
FKAPPA=SQRT(XX) 

DELFN1=.305XKAPPAAXDEFFXFF1XFF11XDELXX3/(JARPRLXFKAPPAXXXSQRT(Z)) 
++(7.0  035E-2)X(2.XKAPPAA+KAPPA)XFF1XDELXX2/(FKAPPAXZ)-.305XFKAPPAX 

+FF1XDELXX3/CJARPRLXZXX1.5) 
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FILE:    VAPBL  FORTRAN      A 

DEL  FN2  =  F*G(  PHI  )*DEL XX3+3.  22XFF1XFKAPPAXDELXX2/C  J  ARPRLXSQRTCZ) ) 
+-.14007*<2.*KAPPAA+KAPPA)*FF1XDEL/FKAPPA 

DELFN3=X/DEFF-FXGl(PHI)XDELXX4/3.-1.61XFKAPPAXFFllXDELxx3/(JARf'RL 
+XSQRT(Z))-.305XKAPPAAXDEFFXFF1XFF12XSQRT(Z)XDELX*3/(JARPRLXFKAFPAX 

+X)+.3  05*KAPPAAXDEFF*DXFNXFF1*FF11*SQRT(Z)*DEL**3/(JARPRLXFKAPPA* 
+XXX2) 

DELFN  =  (DELFN3-DELFN1XDZVDELFN2 

ENDIF 

END 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

FUNCTION  ZFN(PHI,Z,DEL) 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

X 

X      STREAMWISE  DERIVATIVE  OF  VAPOR  MOMENTUM  THICKNESS  PARAMETER 

x 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

REAL  ZFN,PHI,Z,DEL,DEFF,A,K,JARPRL,F 

EXTERNAL  XFN 
EXTERNAL  Fl 
EXTERNAL  Fll 

COMMON/GEOM/A, K, DEFF/PARAM/JARPRL , F 

IFCPHI.EQ.O.D+0)  THEN 

ZFN=0.D+0 
ELSE 

ZFN=( .4<*1X(1.-2.XJARPRLXSQRT(Z)/DEL)-6.XDEFFXF11(PHI)XZ/XFN(PHI)) 
+XXFN(PHI)/F1(PHI)/DEFF 
ENDIF 

END 
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

SUBROUTINE  NUSSELTCNEQ, PHI, Y, YPRIME) 
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

X 

X      STREAMWISE  DERIVATIVE  OF  CONDENSATE  FILM  THICKNESS  USING 

X      A  NUSSELT  ANALYSIS  DOWNSTREAM  OF  VAPOR  SEPARATION  POINT. 

x 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

INTEGER  NEQ 

REAL  PHI,Y(NEQ),YPRIME(NEQ),A,K,DEFF,JARPRL,F 

EXTERNAL  XFN,G,G1 

COMMON/GEOM/A, K, DEFF/PARAM/JARPRL, F 

YPRIME(1)=(XFN(PHI)/DEFF-F/3.XY(1)XX4XG1(PHI))/(FXG(PHI)XY(1)XX3) 

RETURN 
END 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx* 

SUBROUTINE  DRVCFN, PHI, FN1) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

X 

x      SUBROUTINE  TO  NUMERICALLY  DIFFERENTIATE  A  FUNCTION 

K 

*xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
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FILE:  VAPBL     FORTRAN   A 

REAL  FN,PHI,FN1,EPS 
FNl=(FN(PHI+l.E-6)-FN(PHI-l.E-6))/2.E-6 

RETURN 
END 

kxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

SUBROUTINE  TAUFN(PHI,Z, DEL, DZ, DDEL, TAU, DTAU) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
X 

X      SUBROUTINE  TO  DETERMINE  THE  INTERFACIAL  SHEAR  STRESS  AND 

X      ITS  STREAMWISE  DERIVATIVE. 

x 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

REAL  PHI, Z, DEL, DXFN,DZ, DDEL, KAPPA, KAPPA1 , KAPPAA, DKAPPA, 
+DKAPPA1,DKAPPAA,TAU,DTAU,A,K,DEFF,JARPRL,F 

EXTERNAL  XFN,F1,F11 

COMMON/GEOM/A,K,DEFF/PARAM/JARPRL,F 

CALL  DRV(F11,PHI,F12) 
CALL  DRV(XFN,PHI,DXFN) 

KAPPA=DEFF*F11(PHI)XZ/XFN(PHI) 
KAPPA1=JARPRLXSQRT(Z)/DEL 
KAPPAA=  .  0682+.174XKAPPA1 

DKAPPA  =  DEFF/XFNCPHI)X(ZXF12+F11(PHI)XDZ)-DEFF/(XFN(PHI;))XX2X 
+  F1KPHDXZXDXFN 

DKAPPA1=JARPRLX( . 5xDELxDZ-ZXDDEL)/CSQRTCZ)XDEL*x2) 

DKAPPAA= . 174XDKAPPA1 

TAU=6.44XF11(PHI)XSQRT((KAPPAAXKAPPA+KAPPAAXX2)/Z) 
DTAU=6.44X(((KAPPAAXKAPPA+KAPPAAXX2)XF11(PHI)+.5XF1(PHI)X 
+(KAPPAAXDKAPPA+(2.XKAPPAA+KAPPA)XDKAPPAA))/SQRTCZX(KAPPAAX 

+KAPPA+KAPPAAX*2))-.5X(KAPPAAXKAPPA+KAPPAAXX2)*F1CPHI-)XDZ/ 
+SQRTCZXX3XCKAPPAAXKAPPA+KAPPAAXX2) ) ) 

RETURN 
END 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

SUBROUTINE  FCN(NEQ,X, Y, YPRIME) 
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

X 

X     SUBROUTINE  TO  INPUT  DIFFERENTIAL  EQUATION  TO  BE  SOLVED 

x      BY  THE  PROBLEM  SOLVER. 

x 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

INTEGER  NEQ 

REAL  X,Y(NEQ),YPRIME(NEQ) 

EXTERNAL  ZFN,DELFN 

YPRIME(1)=ZFN(X,Y(1),Y(2)) 
YPRIME(2)=DELFN(X,Y(1),Y(2),YPRIME(1)) 

RETURN 
END 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

SUBROUTINE  FCNJ(NEQ,X, Y, DYPDY) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
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FILE:  VAPBL     FORTRAN   A 

*     SUBROUTINE  REQUIRED  BY  PROBLEM  SOLVER,  NO  FUNCTION  IN  ANALYSIS 

INTEGER  NEQ 

REAL  X,Y(NEQ),DYPDYCX) 

RETURN 
END 
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